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FORE WARD 


Automation  of  the  control  of  an  industrial  process  generally 
entails  three  basic  tasks: 

i)  design  of  a  computer  system  (the  computer 
configuration,  the  sensors,  the  transistors, 
and  the  control  software) ; 

ii)  design  of  a  control  strategy  or  control  law 
Hi)  representation,  i.e.,  modeling  of  the  process 

There  is  normally  a  natural  ordering  of  these  tasks;  that  is, 
one  does  not  design  the  computer  system  without  knowledge  of  the  con¬ 
trol  strategy,  and  one  does  not  design  a  control  strategy  without  a 
model  of  the  process. 

Often  the  first  task  is  less  process-dependent  than  the  others, 
and  common  computer  technology  is  applicable.  The  second  task  is 
straightforward  for  certain  simple  models;  for  example,  linear  systems, 
systems  with  restricted  size,  and  systems  with  only  discrete  or  only 
continuous  control  (decision)  variables. 

For  many  complex  industrial  processes,  it  is  the  second  and 
third  tasks  that  limit  realization  of  automated  control.  This  is  true 
for  the  copper  converter  aisle  described  in  this  report.  The  converter 
aisle  process  is  large,  nonlinear,  and  has  mixed  continuous  and  discrete 
control  and  decision  variables.  For  this  type  of  system  work  is  being 
done  on  development  of  representative  models  and  improved  senr.ors,  but 


in 


In  part  this  is 
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work  on  design  of  control  strategies  has  lagged, 
because,  in  concept  at  least,  optimal  control  design  can  be  achieved 
using  standard  techniques  such  as  Bellman's  dynamic  programming. 
Unfortunately,  for  systems  such  as  the  converter  aisle,  standard  con¬ 
trol  design  techniques  lead  to  prohibitively  large  computational 
tasks . 

This  paper  reports  a  new  and  simole  approach  for  computation  of  the 
control  for  a  process  as  typified  by  the  copper  converter  aisle.  The 
resultant  control  is  open  loop  and  locally  optimal.  The  models  used 
to  represent  the  process  and  to  illustrate  resuJts  are  unrealis¬ 
tically  simple,  but  they  serve  to  illustrate  application  of  the 
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ABSTRACT 


The  task  of  dynamic  optimization  consists  of  manipulating  the 
inputs  to  a  dynamic  system  (i.e.,  one  in  which  the  state  varies  with 
time)  so  that  the  system  performs  in  an  advantageous  manner. 

This  paper  presents  a  systematic  technique  for  solving  the  problem 
of  optimally  controlling  a  converter  aisle  in  a  copper  smelter.  The 
converter  aisle  is  distinguished  from  the  usual  dynamic  system  in  that 
some  of  the  control  variables  occur  as  discrete  decisions  while  others 
may  vary  continuously  with  time.  In  this  sense,  the  converter  aisle 
typifies  many  industrial  processes.  The  aisle  is  viewed  as  a  total  sys¬ 
tem  with  the  objective  of  optimizing  overall  performance  as  evaluated 
using  a  mathematical  performance  criterion.  Typical  criteria  reflect 
total  processing  time  and  operating  costs. 

An  essential  step  towards  optimization  is  the  development  of  a  mathe¬ 
matical  model  to  predict  Uie  state  of  the  system.  A  simple  mathematical 
of  the  converter  aisle  is  developed , and  using  this  model,  two  optimiza¬ 
tion  approaches  are  examined  ,  direct  optimization  of  the  total  system 
and  partitioning  the  system  into  interacting  subproblems.  The  parti¬ 
tioned  approach  is  pursued  in  detail  with  techniques  for  solving  the 
optimization  subproblems  presented  and  illustrated  by  numerical  examples. 


IX 
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INTRODUCTION 

The  converter  aisle  in  a  copper  smelter  is  composed  of  a  number 
of  converters  (or  furnaces)  ,  a  supply  of  input  material,  a  receiver  of 
molten  copper,  and  a  crane  (or  cranes)  which  serves  the  converters. 

The  converters  oxidize  the  input  material  to  produce  first  an  intermed¬ 
iate  material,  and  later  molten  copper.  The  crane  delivers  material  to 
the  converters  and  removes  waste  material  and  molten  copper.  Viewed  as 
a  system  to  be  controlled  in  some  advantageous  manner,  the  aisle  has 
an  important  characteristic  in  that  some  of  the  "controls"  are  discrete 
decisions  (for  example,  the  times  that  the  crane  visits  each  converter) 
and  others  are  continuous  variables  (for  example,  oxygen  flow  rate  to 
each  converter) . 

The  application  of  optimal  control  techniques  to  industrial  pro¬ 
cesses  has  commonly  been  limited  to  techniques  for  static  optimization, 
such  as  allocation  and  scheduling  problems  involving  discrete  decisions 
and  "jet-point"  type  control  of  continuous  variables.  Because  of  its 
complexities,  there  has  been  little  attention  given  to  the  copper  aisle. 
Recently  some  papers  have  reported  on  this  process,  but  these  have  dealt 
only  with  parts  of  the  total  optimal  control  problem. 

This  paper  presents  a  technique  for  solving  the  problem  of 
optimal  control  of  the  converter  aisle  viewed  as  a  total  system.  We 
first  establish  a  perspective  of  the  problem  from  a  general  systems  view¬ 
point  and  then  view  optimal  control  of  a  converter  aisle  as  a  general 
systems  optimization  problem.  We  point  out  major  steps  and  approaches 
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to  aisle  optimization.  Two  optimization  approaches  are  presented.  The 
first  is  a  direct  approach  which  results  in  a  problem  that  is  computa¬ 
tionally  large— so  large,  in  fact,  that  a  major  investment  in  time  and 
effort  is  required  to  produce  any  results.  Consequently,  the  direct 
optimization  approach  is  not  used  but  rather,  an  alternative  method 
is  developed  whereby  the  problem  is  partitioned  into  smaller  interacting 
subproblems  which  may  be  solved  with  a  reasonable  effort  in  a  relatively 
short  period  of  time.  The  contribution  of  this  paper  is  the  systematic 
method  of  partitioning  a  large  complex  problem,  such  as  the  converter 
aisle,  into  a  class  of  more  easily  solved  subproblems  and  then  com¬ 
posing  the  solutions  to  the  subproblems  to  form  a  total  optimal  solution. 


CHAPTER  I 


SYSTEMS  CONCEPTS 


Modern  control  theory  offers  a  richly  developed  body  of  techniques 
for  optimally  controlling  a  dynamic  system.  These  techniques  will  be 
increasingly  applied  to  control  of  industrial  processes.  This  section 
establishes  basic  systems  terminology  and  perspective  relating 
to  process  control.  We  note  the  distinction  between  the  physical  pro¬ 
cess  and  the  system  model  and  then  discuss  the  role  of  the  model  in 
design  of  the  optimal  control. 


Associated  with  a  physical  process,  we  identify  certain  (time  vary¬ 
ing)  quantities  as  controlled  variables,  control  variables  (or  decision 
variables),  and  internal  variables.  Controlled  variables  are  those  which 
we  wish  to  observe,  control,  or  regulate  in  some  prescribed  manner,  or 
which  are  used  in  evaluation  of  a  performance  criterion  which  measures 
the  process  performance.  Control  variables  are  variables  to  be  manipula¬ 
ted  or  decisions  to  be  made  so  as  to  achieve  desired  process  performance, 
interna)  variables  are  those  additional  variables  which  are  physically 
related  (directly  or  indirectly)  to  either  the  controlled  or  control  vari¬ 
ables.  As  a  conceptual  abstraction,  the  physical  process  is  viewed  as 
a  system  represented  as  the  block  diagram  in  Figure  1,  where  the  arrows 
represent  flow  of  information.  The  system  input  is  a  vector  of  the  pro¬ 
cess  control  variable,,  and  the  system  output  is  a  vector  of  the  controlled 
variables  of  the  process.  The  state  vector  is  derived  from  the  internal 
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and  controlled  variables  of  the  process;  the  state  variables  are  the 
minimum  set  of  variables  necessary  to  specify  at  any  time  the  value 
of  the  output  given  the  value  of  the  input  at  the  same  time.  The 
state  vector  xs  not  necessarily  unique,  as  it  depends  upon  the  choice 
of  output  variables  and  upon  the  degree  of  approximation  in  specifying 
the  value  of  the  output  vector . 

In  design  of  the  system  control,  it  is  necessary  to  evaluate  the 
output  in  response  to  choices  of  the  input.  Generally  it  is  unfeasi¬ 
ble  or  prohibitively  expensive  to  experiment  with  the  actual  process; 
consequently,  a  model  of  the  system  is  developed.  The  model  consists 
of  two  sets  of  rules  which  specify  the  output  and  the  st  te.  The  output 
equation  is  generally  algebraic  and  it  specifics  the  output  as  a  function 
of  the  state  and  the  input.  Commonly,  the  output  will  be  merely  a  sub¬ 
set  of  the  state  variables.  The  state  equation  depends  upon  the  type  of 
model,  there  being  two  primary  types  of  models,  static  models  in  which 
the  state  does  not  vary  with  time  and  dynamic  models  in  which  the  state 
does  vary  with  time.  In  addition,  a  dynamic  model  is  stationary  if  the 
equations  do  not  vary  with  time.  For  a  static  model,  the  state  equations 
are  generally  algebraic.  Dynamic  model  state  equations  take  a  variety 
of  forms:  for  an  analytic  continuous  model,  i.e.  ,  with  continuous  vari¬ 
ables,  the  state  equations  are  differential  equations  which  specify  the 
time  rate  of  change  of  the  state  variables.  If  the  state  variables  are 
naturally  discrete  or  are  quantized  in  time,  the  discrete  model  may  be 
^  f f srence  equation  or  state  transition  equation  which  specifies  the 
next  state  as  a  function  of  the  current  state  and  input.  Alternately,  a 
state  transition  equation  may  be  given  in  a  tabular  or  graphic  form. 
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Computer  simulation  models  incorporate  analytic,  tabular,  and  partially 
verbal  representationsof  the  model.  Whereas  the  analytic  and  tabula:.- 
models  will  be  directly  useful  in  design  of  optimal  controls,  simulation 
models  are  commonly  used  for  parametric  experimentation  on  the  system. 

The  system  represented  in  Figure  1  is  si, own  as  a  single  block, 
although  systems  are  often  modeled  as  a  set  of  interacting  subsystems, 
for  example,  as  illustrated  in  Figure  2.  System  structure  may  be  natu¬ 
re]  ly  implied  by  identification  of  subsystems  associated  with  the  physi¬ 
cal  process,  or  the  system  structure  may  be  imposed  to  achieve  computa¬ 
tional  advantages.  Decomposition  of  a  system  into  subsystems  has  the  ad¬ 
vantage  that  the  individual  subsystems  are  more  easily  analyzed,  but  on 
the  other  hand,  the  interaction  of  the  subsystems  may  be  difficult  to 
handle.  (it  should  be  kept  in  mind  that  subsystem  input  variables  are 

not  necessarily  physical  quantities,  but  may  represent  decisions,  sched¬ 
ules,  etc.) 

Generation  if  a  system  model  is  a  complex  task  (so  much  so  that 
there  is  danger  of  temporarily  ignoring  the  optimization  objective) .  The 
modeling  task  has  two  basic  steps,  determination  of  the  model  form,  i.e., 
Structure  and  type,  and  data  acquisition  and  numerical  construction  of 
the  model.  The  form  of  the  model  is  usually  motivated  by  physical  knowl¬ 
edge  of  or  familiarity  with  the  process.  Numerical  design  of  the  model 
often  reduces  to  statistical  estimation  of  a  set  of  parameters,  which  in 

turn  involves  subtasks  of  filtering  data  and  estimating  state  variables 
which  are  not  directly  measurable. 

The  role  of  optimal  control  theory  is  to  provide  a  method  for  deter¬ 
mining  dynamic  inputs  which  cause  the  system  to  behave  optimally  i„  the 
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that  a  performance  criterion  J,  which  measures  the  influence  of 
the  inputs  upon  the  system  performance,  will  be  minimized  (or  maxi¬ 
mized)  subject  to  system  constraints.  For  application  of  optimization 
techniques  the  performance  criterion  is  commonly  of  the  form 


J  =  0 

x(t 

! ,  t  | 

1  + 

— '  f ) 

f  | 

<P(X,  U,  T)dT 


(1.1) 


for  a  continuous  dynamic  model 


or 


N 


J  “  I  u,  t] 

t=0 


(1.2) 


for  a  discrete  model,  where  6  and  *  are  functionals  (i.e.,  they  take  on 
scalar  values)  and  tQ  tc  tf  is  the  time  interval  of  interest  or  t  -  0 . 

“  “e  the  diS=rete  timeS  of  although  the  final  time  need  not  be 

fixed.  The  performance  criteria  Equations  1.1  and  1.2  are  functions  of 

the  system  trajectory  U(t>  ,  „|t)},  that  is,  the  time  history  of  the 
state  and  associated  control.  The  system  constraints  reflect  the  physical 
or  economic  constraints  on  the  process.  Given  the  system  model,  perform¬ 
ance  criterion,  and  constraints,  th e  optimisation  problem  may  be  represen¬ 
ted  by  the  diagram  in  Figure  3,  whrch  sb  aid  be  distinguished  from  the 
system  block  diagram.  (The  output  of  the  problem  diagram  is  the  optimal 
solution  u°  which  is  the  input  to  the  system.) 

Depending  upon  the  model,  the  criterion,  and  the  constraints,  there 
are  a  number  of  techniques  available  for  determining  u«.  Commonly,  u° 
is  found  as  a  control  policy  u°  -  „"(t)  so  that  the  control  is  "open 
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loop.  In  contrast,  for  practical  implementation,  it  is  desired  that 
the  control  be  a  feedback  lav/  u°  =  u°  [x,t],  a  function  of  the  measured 
(cr  estimated)  state,  as  illustrated  in  Figure  4.  In  some  cases,  the 
feedback  law  may  be  directly  determined,  but  more  often  a  neighbor 
optimal  controller  is  used.  In  this  case,  the  open  loop  u°(t)  and  cor¬ 
responding  x° (t)  are  calculated  and  a  simple  controller  is  designed  to 
keep  the  system  "near"  the  desired  state  in  spite  of  random  disturbances 
which  cause  deviation  from  desired  performance.  We  now  apply  these 
systems  concepts  to  the  copper  aisle  process. 


CHAPTER  II 


CONVERTER  AISLE  SYSTEM 


«a  will  now  formulate  the  converter  aisle  as  a  general  system, 
in  order  to  identify  the  input  variables,  state  variables  and  output 
variables,  a  sequence  of  events  in  a  converter  aisle  is  briefly  de¬ 
scribed.  For  the  purposes  of  illustration,  let  an  example  converter 
aisle  consist  of  one  reverberatory  furnace,  three  converter  furnaces 
and  a  holding  furnace  as  shown  in  Figure  5. 

Liquid  matte  (FeS  and  C^S)  is  poured  into  the  mouth  of  a  con¬ 
verter  from  a  ladle  which  is  manipulated  by  an  overhead  crane,  often 
cold  material,  including  matte  shells,  scrap  copper,  etc,  is  added  with 
the  initial  charge  of  matte.  The  converter  is  then  rolled  into  the 
blaming  position,  submerging  the  air  inlet  openings  (tuyeres)  below  the 
surface  of  the  molten  hath  and  air  is  blown  into  the  bath,  thus  remov¬ 
ing  by  oxidation,  or  slagging,  the  undesirable  elements  <Fe,s)  which 
are  in  the  bath.  Flu: r  is  added  through  an  overhead  hopper  or  by  means 
Of  a  conveyor  belt  which  dumps  the  flux  into  the  hood  of  the  converter. 
Blowing  is  continued  long  enough  to  use  up  the  flux  and  to  form  slag 
which  is  periodically  skimmed  (poured)  from  the  converter  by  turning  it 
out  of  the  blowing  position.  More  matte,  flux  and  cold  material  are 
added  and  then  the  converter  is  turned  back  to  the  blowing  position 
and  blowing  is  resumed.  These  partial  blows  continue  until  all  the  iron 
sulfide  has  been  oxidised  and  the  converter  contains  essentially  pure 
cuprous  sulfide  (Cu2s) .  Several  of  these  matte  charges  must  be  slagged 
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before  there  is  enough  Cu2S  accumulated  to  proceed  with  the  finish  or 
copper  blow.  As  soon  as  all  of  the  sulfur  in  the  Cu2S  has  been  oxidized, 
the  metal,  blister  copper3  is  poured  out  and  the  converter  is  available 
to  begin  another  charge  sequence. 

As  the  converters  process  material,  the  crane  moves  along  the 
aisle  track  to  deliver  and  remove  material  in  response  to  needs  of  the 
converters.  The  crane  obtains  matte  and  cold  material  at  the  reverber¬ 
atory  furnace  site  and  delivers  slag  to  the  reverberatory  furnace  as 
it  delivers  blister  copper  to  the  holding  furnace,  which  is  assumed  to 
be  located  near  the  reverberatory  furnace.  Generally,  each  service  to 
a  converter  requires  several  trips  by  the  crane.  Also,  when  a  con¬ 
verter  is  being  serviced,  it  does  not  blow  between  visits  by  the  crane. 

The  above  descriptions  suggest  two  activities  that  characterize 
the  converter  aisle:  (1)  material  handling  by  the  cranes  and  (2)  charge 
processing  by  the  converter  furnaces.  The  behavior  of  the  crane  is 
characterized  by  its  location  and  the  full  or  empty  status  of  the  ladle 
which  we  will  establish  together  as  a  two-dimensional  state  vector.  The 
input  variable  which  affects  the  state  of  the  crane  is  a  decision  which 
specifies  the  position  to  which  the  crane  will  move.  The  result,  or  out¬ 
put,  of  the  crane  activity  is  the  transfer  of  material  into  or  out  of  the 
converters.  From  this  it  is  apparent  that  discrete  material  inputs  to 
the  converters  are  directly  associated  with  the  activity  of  the  crane. 

The  second  activity  which  characterizes  the  behavior  of  the  con¬ 
verter  aisle  is  the  processing  of  matte  by  the  converter  furnaces .  The 
input  quantities  that  affect  the  state  of  the  converter  can  be  measured 
in  a  straightforward  way.  Air  flow  into  and  exhaust  gases  from  each 
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converter  can  be  measured  continuously  and  the  material  input  and  out¬ 
put  of  the  converter  can  be  weighed  and  the  composition  analyzed 
with  some  degree  of  confidence.  For  the  slag  phase,  essential  inputs 
indicated  by  the  previously  described  events  are: 

u^  =  Matte  weight 
u^  =  Flux  weight 
u^  =  Cold  material  weight 
u^  =  Oxygen  input  rate 

In  contrast  to  the  input  variables,  the  activities  directly  asso¬ 
ciated  with  the  progress  of  each  converter  are  difficult  to  measure  dir¬ 
ectly.  The  state  variables  chosen  to  characterize  the  converting 
process  are  somewhat  arbitrary  and  certainly  do  not  describe  every  ac¬ 
tivity  of  the  process  exhaustively.  The  state  variables  are  chosen  in 
accordance  with  the  information  needed  to  adequately  and  realistically 
describe  the  observable  and  measurable  activities  of  the  process.  Since 
the  function  of  the  converter  is  to  remove  undesirable  elements  from  the 
copper  bearing  matte,  the  relative  amount  of  FeS  contained  in  the  bath 
at  any  time  is  of  interest  during  the  slag  phase.  Similarly,  the  amount 
of  Cu^S  provides  information  concerning  the  progress  of  the  finish  phase. 
The  amount  of  slag  produced  would  have  an  effect  on  the  requirement  for 
skimming.  The  copper  accumulated  is  of  fundamental  importance  since  that 
is  the  objective  of  the  converter  furnace.  The  bath  temperature  would 
surely  have  an  effect  on  the  reaction  rates  within  the  process.  Sum¬ 
marizing,  a  typical  set  of  state  variables  required  to  describe  the  pro<- 
cess  activities  during  the  slag  phase  are: 
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=  FeS  concentration 

*2  =  ('U2^  weight 
x^  =  Flux  weight 
x^  =  Slag  weight 
x  =  Temperature 

D 

x6  =  Magnetic  (Fe  0  )  weight 


The  finish  phase  is  similar  in  concept  to  the  slag  phase,  except 
that  during  the  finish  phase,  the  bath  consists  essentially  of  Cu2S  which 
is  being  converted  to  blister  copper. 

The  discrete  decision  variables  of  the  system  are: 

!•  Crane  movement  decision 

2.  Converter  schedule 

3.  Batch  material  inputs 

Of  course,  these  are  not  independent  decisions.  The  continuous  control 
variables  are  the  continuous  inputs  such  as  oxygen  blow  rate  and  possi¬ 
bly  flux  addition  rate.  The  rules  for  determining  the  time  history  of 
the  state  variables  are  discussed  in  detail  in  the  succeeding  chapter 
on  modeling.  The  copper  production  rate  is  specified  as  an  output  vari¬ 
able  as  a  matter  of  convenience  in  the  development  of  the  optimization 
objective . 

One  objective  of  optimization  of  the  converter  aisle  system  is  to 
specify  the  control  values  which  will  influence  the  system  to  achieve 

the  greatest  possible  value  for  the  output,  where  the  output  performance 
criterion  is  expressed  by: 


J  = 


3 

I 

i=l 


Cu. 

l 


(2.1) 
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where 

Cui  =  Blister  copper  produced  per  charge  cycle 
in  converter  i 

Ti  =  Oration  of  a  charge  cycle  for  converter  i 

The  charge  cycle  time  is  taken  to  be  the  time  to  complete  a  charge  cycle 
on  a  converter  as  shown  in  Figure  6.  Crane  movement  decisions  influence 
the  charge  cycle  time  particularly  when  two  or  more  converters  are  in 
conflict  for  service.  The  converter  furnace  schedule  affects  the  charge 
cycle  time  by  adjusting  the  time  staggering  between  converters.  The 
batch  material  inputs  influence  the  charge  cycle  times  and  the  blister 
copper  production,  while  continuous  material  inputs  affect  the  blowing 
times.  Thus,  we  observe  that  the  performance  criterion  of  Equation  2.1 
is  an  implicit  function  of  the  state  and  control  variables. 

We  next  turn  to  establishing  a  system  model. 
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Figure  5  EXAMPLE  CHARGE  CYCLE 
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CHAPTER  III 

CONVERTER  AISLE  MODELS 

Motivated  by  the  ideas  presented  in  the  previous  chapter,  we 
will  develop  a  model  of  each  subsystem  of  Figure  7.  The  emphasis  of 
this  section  is  on  modeling  concepts  rather  than  the  exhaustive  deri¬ 
vation  of  a  specific  model.  Thus,  the  models  that  are  developed  are 
purposely  simple  with  the  assumptions  used  clearly  stated. 

It  is  not  possible  to  specify  general  rules  by  which  system 
models  are  formulated,  but  a  number  of  guiding  principles  can  be  stated: 

1.  The  system  will  be  organized  in  subsystems  so  as  to  sim¬ 
plify  the  specification  of  interactions  within  the  system. 

The  subsystems  will  be  characterized  by  one  type  of  variable 
(i.e.,  discrete  or  continuous). 

2.  The  model  will  include  only  those  variables  which  are  rele¬ 
vant  to  the  optimization  objective. 

3.  The  accuracy  of  the  model  will  be  limited  by  the  data  used 
to  develop  the  model . 

4.  The  model  must  be  capable  of  reproducing  the  behavior  of 
the  system  within  an  acceptable  accuracy. 

Crane  System  Model 

The  physical  structure  of  the  converter  aisle  suggests  a  crane 
routing  model.  The  crane  state  will  be  a  two-vector  P_.  which  represents 
the  state  after  the  jth  crane  move: 

r*, . 

D 

j 


P  . 
“3 


(3.1) 
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The  first  component  ?1  represents  the  crane  location,  which  can  have 
four  possible  values 

=  0  crane  at  reverberatory  furnace  or 
holding  furnace 

P^  =  1  crane  at  converter  1 

P^  =  2  crane  at  converter  2 

P^  =  3  crane  at  converter  3 

The  second  component  P2  represents  the  ladle  condition  during  the 

move,  where 


P2  =  0  ladle  empty 
P2  =  1  ladle  loaded 


The  crane  is  modeled  by  a  state  transition  equation  represented  as 


-  !l6v  di) 


(3.2) 


where  P  is  the  current  state,  d  is  the  crane  move  decision,  and  P 

J  -j+1 

is  the  next  state.  P.j+j_  is  determined  by  d^  subject  to  the  following 
assumed  constraints: 


Only  one  ladle  of  material  will  be  input  or  output  from  a 
converter  at  each  visit. 


2.  Only  an  empty  ladle  will  be  moved  to  a  converter  that  needs 
to  be  skimmed  or  emptied. 

3.  Only  a  loaded  ladle  will  be  moved  from  a  converter  that  has 
been  skimmed  or  (partially)  emptied. 

4.  Only  a  loaded  ladle  will  be  moved  to  a  converter  which  re¬ 
quires  input  material,  provided  it  does  need  output  service 

5.  Only  an  empty  ladle  will  be  moved  from  a  converter  after 
delivery  of  input  material. 
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With  these  assumptions,  it  is  seen  that  choice  of  the  value  of  the 
decision  d.  requires  knowledge  of  the  converter  states;  and  given 
that  knowledge,  the  decision  d.  is  a  scaler  value,  i.e.,  it  is  neces¬ 
sary  only  to  choose  (P^  subject  to  the  constraints  and  then  (P  ) 

2  j+1 

is  uniquely  determined. 

For  a  given  crane  decision,  the  additional  information  required  to 
model  the  crane  is  the  time  required  to  move  from  the  current  location 


to  the  next  location.  This  transition  time  information  will  be  given  in 
the  form  of  Tables  1  and  2,  where  it  is  assumed  that  the  transition 
time  includes  the  loading  time  at  the  reverberatory  furnace  and  the  load¬ 
ing  or  unloading  time  at  the  converter  to  which  the  crane  moves.  The 
blanks  in  the  tables  indicate  illegal  moves. 

As  noted  above,  crane  movements  are  made  in  response  to  material 
input/output  needs  of  the  converters,  and  the  moves  in  turn  effect  the 
states  of  the  converters.  The  constraints  of  converter  requirements  on 
the  crane  decision  can  be  described  in  terms  of  the  converter  state 
variables,  and  the  effect  of  a  crane  move  on  the  converters  can  be 
modeled  by  including  the  crane  variables  in  the  converter  models.  (This 
is  done  in  Chapter  V  in  the  direct  optimization  approach.)  However, 
even  if  the  converters  are  not  explicitly  modeled,  it  may  still  be  desir¬ 
able  to  deal  with  the  crane  subsystem.  For  this  case,  to  model  the 
crane  interaction  with  the  converters,  we  model  requirement  variables  as 
an  input  queue  and  an  output  queue.  The  input  queue  is  a  vector 


—in  = 


m. 


m„ 


m. 
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TABLE  1 

STATE  TRANSITION  TIMES,  (p =  1 

V  2Vl 

State 


TABLE  2 

STATE  TRANSITION  TIMES,  fp  ^  0 

V  2/j+l 
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where  m.  is  the  number  of  ladles  to  be  input  to  the  ith  converter. 


The 


output  queue  is  expressed  by 

^out 


Si  is  the  number  of  ladles  of  material  to  be  removed  from  the  ith 
converter.  Clearly,  the  queue  variables  can  have  only  non-negative  in¬ 
teger  values,  in  terms  of  the  queues,  the  crane  moves  only  in  response 
to  nonzero  queue  values  or  to  return  to  the  reverberatory  or  holding 
furnace.  Each  time  the  crane  moves  to  a  converter,  the  appropriate 
queue  variable  for  that  converter  decreases  by  one.  Each  time  the  crane 
moves  from  a  converter,  P2  must  change,  i.e.,  (P2) j+l  +  (p^  ..  A  crane 

move  to  the  reverberatory  furnace  or  holding  furnace  does  not  change  the 
queues . 

In  addition  to  the  changes  in  queue  variables,  the  state  of  the 
converter  must  change  since  material  is  being  removed  or  added.  The 

Change  in  COnVerter  state  —  s ervi ce  is  given  by  the  following 

relationship 

x  (k+1)  =  x  (k)  +  e  m . 

1  1  1 

x  (k+1)  =  x  (k)  +  e  m. 

*  ^  2  i 

x  (k+1)  =  x  (k)  +  e  e  m. 

J  J  3  2  i 

x^ (k+1)  =  0 

x5(k+D  =  (x5(k)x2(k)  +  (e^.  +  e2m.  +  e^m.)  e^./ 
[x^k+1)  +  x2(k+l)  +  x3  (k+1)  ] 
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x  (k+1)  =  0 
6 

where  the  subscripts  (k+1)  denote  the  time  following  crane  service 
and  (k)  the  time  preceeding  crane  service-  The  parameters  e  ,  e  f 
e^,  e^  are  defined  as 

e1  =  weight  of  FeS  per  ladle  of  input 

e2  =  weight  of  Cu2S  per  ladle  of  input 

e3  =  numt>er  of  tons  of  flux  required  per  ton  of  FeS 

e4  =  average  temperature  of  input  material 

The  model  of  temperature  change  is  a  simple  linear  extrapolation  and 
assumes  there  is  no  change  in  bath  temperature  due  to  radiation  losses 
during  crane  service. 

The  above  crane  service  state  transition  equations  for  the  ith  con¬ 
verter  may  be  summarized  by 


x(k+l)  =  x(k)  +  A  xfejirt)  (3.3) 

where  e  is  a  vector  of  model  parameters  defined  above  and  m  is  the  num¬ 
ber  of  ladles  of  input  material. 


The  number  of  ladles  of  material  to  be  removed  is  a  function  of  the 
amount  of  material  to  be  removed  and  the  ladle  capacity.  For  the  ith 
converter  this  may  be  expressed  by 


n . 

l 


x4(k) 

W 


for  the  slag  and  by 
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n . 
1 


(-.8)x2(k) 

w 


+ 


for  the  finish  phase,  where: 


x4  (k)  =  amount  of  slag  at  end  of  slag  blow 


x  (k)  =  amount  of  Cu  S  at  the  start  of  the  finish 
blow  2 


and  where 


w  -  weight  capacity  per  ladle  (assume  slag  and 

copper  are  of  equal  density) 

+ 

t  •  ]  denotes  the  nearest  equal  or  greater  integer. 


The  crane  service  time  between  any  two  consecutive  blows  is,  of 
course,  a  function  of  the  routing  used  in  performing  service.  If  a 
converter  is  serviced  exclusively  (i.e.,  the  crane  is  not  shared  with 
another  converter)  then  the  service  time  is  given  by 


Ts(i'j)  =  (tio  +  +  trj 


oi  i  '  oi 
.th 


10 


m. 


(3.4) 


where  TS(i,j)  is  the  time  to  service  the  jt  blow  on  the  i*^  converter. 

This  modeling  of  the  crane-converter  interaction  completes  a  model 
of  the  crane  subsvstem. 


Converter  Furnace  Model 

We  recognize  from  Chapter  II  that  the  activities  which  character¬ 
ize  the  behavior  of  the  converter  furnace  involve  two  major  reactions. 
The  slag  phase  is  characterized  by  the  elimination  of  iron  sulfide  from 
the  matte  charge  as  governed  by  the  following  chemical  equations: 

2FeS  +  30 2  -+  2FeO  +  2S02  AH  =  -105,000  Cal/mole 

2Fe0  +  S.O  -*  2Fe0  •  S.0„  AH  =  -306 

12  i  2 


Cal/mole  (3.5) 
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6FeO  +  °2  2Fe3°4  AH  =  “133/000  Cal/mole 

3Fe304+FeS  ■+  lOFeO  +S02  AH  =  +110,000  Cal/mole 

where  AH  is  the  heat  of  reaction  at  molten  temperatures.  The  second 
phase  is  described  by  the  reaction 

Cu2S  +  °2  "*■  2Cu  +  S02  AH  =  _56'000  Cal/mole  (3.6) 

In  practice,  these  equations  have  been  used  for  many  years  to  estab¬ 
lish  material  balance  and  heat  balance  relationships.  These  equations 
do  not,  however,  offer  any  information  regarding  the  behavior  of  the 
process  in  the  course  of  time.  It  is  an  interest  in  controlling  the 
processing  time  of  the  converter  furnace  that  motivates  us  to  develop 
a  dynamic  model  which  will  predict  the  value  of  the  state  variables  at 
any  time.  Defining  the  slag  phase  state  variables  as  in  Chapter  II, 
we  wish  to  formulate  a  slag  phase  model  using  oxygen  input  rate  as  the 
only  continuous  control . 

During  each  blow,  the  process  involves  only  continuous  vari- 
bles ,  so  that  it  may  be  modeled  as  a  vector  differential  equation  of 
the  form 

dx 

—  =  x  =  f(x,  u)  (3.7) 

where  x  is  the  state  vector  defined  in  Chapter  II  and  u  is  the  oxygen 
input  rate.  The  function  f_  is  in  general  a  nonlinear  function  of  the 
arguments  x,  u  and  it  is  assumed  to  be  stationary.  Because  of  limited 
theoretical  knowledge  about  the  dynamics  cf  copper  converter  furnaces, 
we  will  approximate  the  function  _f  with  an  empirical,  stationary, 
linear  function.  We  offer  no  experimental  justification  for  such  a 
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simple  model  but  use  it  to  illustrate  a  modeling  technique  and  to  serve 
as  a  simple  model  in  optimization  examples.  While  the  assumption  of 
linearity  may  seem  unjustified,  it  does  have  practical  merit.  A  non¬ 
linear  system  such  as  Equation  3.7,  operating  "near"  a  known  trajec¬ 
tory  (xo(t) ,  UQ(t)}  may  be  linearized  by  the  expansion. 


,  3f , 
f(x  ,  u  )  +  _= 
—  — o  o 

8x 


(x  -  X  )  + 
—  — o 


X  ,  u 
-o  o 


3f 

3u 


(u  -  u  ) 
o 

u 

o 


Evaluating  the  partial  derivatives  as  indicated  and  assuming  equality 
yields  an  equation  of  the  form 


x 


A(t)x  +  B(t)u  +  C(t) 


This  equation  represents  a  nonstationary  linear  system.  A  further 
approximation  to  this  system  is  made  by  assuming  the  system  is  station¬ 
ary  over  an  interval  of  time.  For  the  slag  phase  model,  we  assume  the 
system  is  stationary  over  all  of  the  blows,  thus  yielding  an  approxi¬ 
mate  model  with  form 


>c=Ax  +  Bu+C 

where  A,  B,  and  C  are  constant  matrices  which  will  be  evaluated. 

To  further  simplify  the  model,  we  assume  that  x  (Cu  S  weight) 
is  cons  cant  during  slag  blows  and  that  all  other  bath  constituents 
vary  in  proportion  to  the  FeS  oxidation  rate.  (These  assumptions  are 
supported  by  observation  of  converter  data.)  With  these  assumptions, 
the  slag  blow  state  equations  become 
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X1  '  A11X.  +  A12X2  + 


+  A16X6  +  V  +  C1 


x2  ’  0 
X3  ‘  Vl 


X.  =  K.X, 

4  4  1 


(3.9) 


X5  "  A51X1  +  A52X1  * 


*  A56X6  +  V  +  C5 


X6  '  Vl 


The  unknown  constants  in  Equation  3.9  will  be  determined  empira- 
cally  based  upon  observation  and  estimation  of  converter  state  trajec¬ 
tories.  Direct  measurements  of  many  of  the  state  variables  associated 
with  the  converter  furnace  are  difficult  to  obtain.  Consequently,  some 
of  the  data  is  obtained  indirectly,  such  as  through  the  use  of  material 
balance  relationships.  For  example,  smoothed  trajectories  of  the  state 
variables  are  shown  in  Figure  8  for  a  typical  converter  furnace  with 
fixed  oxygen  input  rate  during  each  blow.  The  original  data  was  obtained 
by  direct  measurement  of  S02  discharge  rate,  oxygen  input  rate,  tempera- 
ture ,  and  initial  material  inputs  and  by  computation  of  the  trajectories 
of  the  concentration  of  FeS  and  weights  of  flux,  slag,  and  Fe.^ 
using  the  reaction  Equations  3.5.  Denote  the  data  points  in  Figure  8 

by  x  (t,  ,)  and  the  slope  at  each  data  point  as  x  (t.  ,)  ,  where  k  denotes 
K  1J  k  i] 

the  kth  state  variable  and  t _  the  ith  discrete  time  in  the  jth  slag  blow. 

To  evaluate  the  constants,  in  particular  the  B  coefficients,  we 

k 

require  a  collection  of  data  for  different  values  of  u.  However,  using 
only  the  data  in  Figure  8,  a  very  simple  initial  model  can  be  obtained 
by  arbitrarily  choosing 


WEIGHT  OF  FLUX  (TONS)  WEIGHT  OF  Cu2S  (TONS)  CONCENTRATION  OF  FeS  (PERCENT) 
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TIME  (MINUTES) 


Figure  8  EXAMPLE  STATE  VARIABLES  (SLAG  PHASE) 


Figure  8  (cent)  EXAMPLE  STATE  VARIABLES  (SLAG  PHASE) 
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_  x,  (t,  .) 

Bk  -  Ji.-1;  '  k  ~  1 '  5  (3.10) 

u(t.  ,) 

13 

where  the  superscript  bar  denotes  the  average  over  all  sample  times 
ti;..  Equation  3.10  has  some  physical  justification  in  that  oxygen  input 
rate  is  dominant  in  determining  FeS  oxidation  rate.  To  evaluate  the 
remaining  constants,  arrange  the  data  in  the  matrix  forms 


xl(tll> 

•  •  • 

x6(tn> 

! 

1  | 

w 

X2(t21> 

•  •  • 

X6,t21> 

1 

XI(t43J 

X2(t43) 

•  •  • 

X6<t43l 

1  i 

and 


4  * 


xk(tn>  -  V(tu> 
W  -  V(t21> 


\(t43>  ‘  V'V 


,  k  =  1,  5 


For  a  least  sum-square-error  fit  of  the  data, 
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-1 

=  £  £  XfcY^  /  k  =  1,5 


(3.11) 


where  superscript  t  denotes  transpose.  The  constants  K  ,  K  ,  and  K 

3  4  6 

in  Equation  3.9  may  be  estimated  as 


h  °  W  -  «  =  3,  4,  6 


(3.12) 


To  complete  the  model,  the  system  constraints  are  simply 


0  <  u(t)  <.  u 
—  —  max 


0  1  xk ( 1  xk  max  /  k  1 , . .  .  ,  6 


(3.13) 


where  and  Xk  max  are  fixed  constants. 

Together,  Equations  3.9-3.13  constitute  a  simple  initial  slag 
blow  moael.  The  model  has  validity  only  during  the  blows  and  only  for 
state  and  control  values  in  the  range  of  data  used  to  derive  the  model 
coefficients.  Between  slag  blows,  oxidation  stops  and  the  converter 
state  changes  in  accordance  with  equation  2.3  as  previously  discussed. 
Thus,  between  blows,  the  crane  actions  determine  the  change  in  the  con¬ 
verter  state.  When  a  sufficient  amount  of  Cu2S  has  been  accumulated 
and  when  the  FeS  has  been  removed,  the  finish  phase  is  started. 

A  finish  phase  model  can  be  developed  in  manner  similar  to  that 
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used  for  the  slag  blow  model,  but  with  appropriate  changes  in  the  state 
variables.  A  typical  state  trajectory  for  the  finish  phase  is  shown 
in  Figure  9.  The  temperature  during  the  finish  phase  remains  relatively 
constant  for  a  given  oxygen  input  rate.  As  expressed  by  Equation  3.6, 
there  is  only  one  principle  reaction  involved,  thus  the  state  vector  will 
be  approximated  by  a  scalar  as 


X2  =  CU2S  weight  (3.14) 

As  in  the  slag  phase  model,  the  finish  phase  model  may  be  expressed  as 
a  differential  equation  of  the  form  of  Equation  3.7. 


d'k2  =  f  ( x  ,  u) 
dt 


(3.15) 


where  f  is  to  be  determined.  The  curve  of  Figure  9  suggests  a  non. 
linear  exponential  form  for  the  function  f 


f (x2,u) 


2  a  ute 


a  ut 


(3.16) 


where  a  is  a  constant  to  be  determined  and  t  is  the  independent  variable 
time.  The  choice  of  an  exponential  function  is  based  simply  on  the 
fact  that  data  best  fit  this  form  of  curve.  Clearly,  the  choice  of  f 
is  not  unique  but  represents  a  judgment  based  on  observation  and  experi¬ 
ence.  Combining  Equations  3.15  and  3.16,  we  have  the  finish  phase  model 
given  by 

x2  =  -2  a  ute  (3.17) 

the  constant  a  can  be  determined  by  solving  equation  3.15  and  fitting 
the  solution  to  the  data  in  Figure  9  using  a  standard  curve  fitting  tech- 


MINUTES 
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nique,  similar  to  the  least  squares  method  used  in  modeling  the 
slag  phase. 

The  model  illustrated  in  this  section  is,  of  course,  quite 
simple,  and  much  remains  to  be  done  in  modeling  of  copper  conversion 
dynamics.  Towards  this  objective,  more  complex  models  have  been  con¬ 
sidered,  for  example,  Nenonen  and  Paqurek  [8]  and  Niemi  and  Koskineu 
[9].  However,  this  simple  model  will  provide  a  basis  for  the  optimi¬ 
zation  concepts  which  follow. 
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CHAPTER  IV 

SIMULATION  STUDIES 

In  this  chapter  a  simulation  of  the  simple,  single  crane,  three 
converter  aisle  shown  in  Figure  5  is  described.  The  purpose  of  the 
simulation  is  twofold;  first,  to  validate  the  model,  and  second,  to 
provide  the  motivation  for  studying  the  converter  aisle  as  an  optimiza¬ 
tion  problem.  In  the  simulation  study  opportunities  for  optimal  control 
are  identified  and  a  measure  of  possible  improvement  in  performance  is 
established.  The  simulation  program  is  briefly  described,  some  results 
are  presented,  and  a  summary  of  promising  optimization  opportunities  is  given 

Simulation  Program 

This  simulation  study  is  conducted  with  the  use  of  an  IBM  1800 
Process  Control  Digital  Computer  and  a  FORTRAN  based  simulation  prog¬ 
ram  called  GASP  II  (General  Activity  Simulation  Program)  [15].  Two 
principle  activities  of  the  converter  aisle  system  are  identified  as 
crane  service  and  processing  of  material  by  the  converters  ^blowing) . 

The  crane  service  activity  is  performed  by  an  entity  of  the  system, 
namely  the  crane,  while  the  blowing  activity  is  performed  by  another 
entity,  the  converter  furnace.  Each  of  these  entities  have  certain 
attributes  which  characterize  their  behavior  as  shown  in  Table  3. 

Clearly  a  mixture  of  continuous  and  discrete  event  models  is 
found  at  the  boundary  between  the  crane  subsystem  and  the  converter  fur¬ 
nace  subsystem.  Two  types  of  events  are  specified.  A  type  1  event  is 
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defined  at  the  endpoint  of  a  trajectory  which  satisfies  the  differ¬ 
ential  equations  of  the  converter  furnace  model  (i.e.,  event  type  1 
occurs  at  the  end  of  a  blow) .  This  event  represents  a  discrete  deci¬ 
sion  point.  The  type  2  event  is  defined  at  the  completion  of  crane  ser¬ 
vice.  This  event  represents  a  continuous  decision  point  since  it  is 
necessary  at  this  event  time  to  decide  oxygen  input  rate  for  the  blow 

that  is  to  follow.  Figure  10  illustrates  the  event  model  used  in  the 
simulation. 

Since  there  are  three  converters  operating  simultaneously,  the 
program  maintains  an  event  file  as  shown  in  Table  4.  The  time  a  partic¬ 
ular  event  will  occur  is  stored  in  the  event  time.  On  each  iteration 
in  the  program,  the  smallest  event  time  is  determined  and  the  event 
code,  converter  number  and  blow  number  are  removed  from  the  event  file. 
At  this  point,  another  file  is  accessed  to  determine  the  duration  of  the 
activity  associated  with  the  event  code.  The  file  containing  the  activ¬ 
ity  duration  is  the  activity  file  as  illustrated  in  Table  5.  The  pro¬ 
gram  (EVNTS)  which  stores  and  retrieves  information  from  these  files  is 
shown  in  Appendix  1. 

The  event  code  determines  which  one  of  two  model  subprograms  is  to 
be  executed.  It  the  event  code. is  1,  then  the  event  marks  an  end  of 
blow,  and  crane  service  is  called  for.  The  crane  service  subprogram 
(CRNSV;  appendix  1)  calculates  the  change  in  state  due  to  the  action  of 
the  crane  and  also  the  duration  of  the  crane  service  activity,  it  is 
assumed  that  once  crane  service  has  been  started  on  a  converter,  it  will 
continue  until  the  needs  of  that  converter  have  been  satisfied  (i.e. 
time-sharing  of  the  crane  is  not  allowed) . 


Thus  the  crane  service  mod- 
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State  Prior  to  Crane  Servi 


Discrete 
Material 
Input  — 


EVENT  TYPE  2 

X(k  +  1) 


►  Crane 
Service 
T  i  me 


State  Prior  to  Glowing 


Continuous 
Oxygen 
Input  - — 


CONVERTER 

FURNACE 

MODEL 


EVENT  TYPF  1 

x(tf) 


♦Blow 

Duration 

Cto5tf] 


Figure  io  SYSTEM  EVENT  MODEL 
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TABLE  4 

TYPICAL  EVENT  FILE 


EVENT 

TIME 

EVENT 

CODE 

CONVERTER 

NUMBER 

BLOW 

NUMBER 

0.0 

1 

1 

3 

50.0 

2 

2 

4 

75.0 

1 

_ 

3 

1 

TABLE  5 

TYPICAL  ACTIVITY  FILE 


ACTIVITY 

CONVERTER 

DURATION 

NUMBER 

70.0 

23.0 

2 

50.1 

_ , 

3 
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el  is  given  by  Equations  3.3  and  3.4.  The  time  for  the  next  event  on 
the  particular  converter  is  calculated  by  adding  the  activity  dura¬ 
tion  to  the  current  time.  The  event  code  is  changed  to  an  end  of  crane 
service  and  a  blow  is  initiated.  The  converter  furnace  subprogram 
( CONVT ;  appendix  1)  calculates  the  change  in  state  due  to  the  activity 
of  the  converting  process.  The  converter  furnace  model  as  given  by' 
Equations  3.9  for  the  slag  activity  duration  is  the  time  required  for 
the  converter  to  transfer  the  state  from  the  initial  condition  at  the 
end  of  crane  service  to  a  specified  terminal  state  constrained  by 
W  1  0  for  i  =  1,  2,  3,  4,  6 

x_(tj  <  2250 °F 
of  — 

Violation  of  any  of  the  above  terminal  constraints  will  terminate  the 
blow  simulation.  In  addition,  the  last  slag  blow  and  finish  blow  have 
special  terminal  conditions.  At  the  end  of  the  final  slag  blow  X;L(t  )  = 
0  and  at  the  end  of  the  finish  blow  x2(tf)/x2(0)  is  specified. 

The  model  parameters  and  assumed  oxygen  input  rates  are  found  in 
Appendix  2. 


Simulation  Results 

Figure  11  shows  some  sample  state  trajectories  obtained  from  a  sol¬ 
ution  of  the  slag  phase  model  and  a  comparison  of  interval  1  trajector¬ 
ies  from  the  model  with  interval  1  trajectories  from  observed  data.  Fig¬ 
ure  12  illustrates  a  comparison  of  the  state  trajectory  obtained  from 
the  solution  of  the  finish  phase  model  and  the  state  trajectory  from 
actual  observation.  The  simulation  results  and  actual  results  are  nc*- 
directly  comparable  because  some  of  the  parameters  associated  with  the 
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Figure  ii  STATE  VARIABLES  (SLAG  PHASE)  COMPUTED  FROM  MODEL 
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TIME  MINUTES 


Figure  li  (cont)  STATE  VARIABLES  (SLAG  PHASE)  COMPUTED  FROM  MODEL 
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actual  data  are  unknown  and  thus  are  assumed  for  the  model .  These 
assumptions  lead  to  some  differences  in  initial  conditions  between 
the  actual  data  and  model  results  as  observed  in  Figure  ,lx.  However, 
comparisons  made  between  actual  observation  and  model  results  showed 
considerable  similarity;  and,  therefore ,  this  model  will  serve  as  a  nomi¬ 
nal  starting  point  for  further  optimization  studies. 

The  results  of  six  simulation  runs,  each  with  a  different  assump¬ 
tion  or  initial  condition  show  a  distinct  difference  in  process  perform¬ 
ance.  The  model  parameter  values  for  each  run  are  shown  in  Append¬ 
ix  2.  The  first  two  simulation  runs  are  designed  to  illustrate  a  dif¬ 
ference  in  process  performance  due  to  two  different  material  input  deci¬ 
sions  as  given  in  Table  6.  The  results  of  this  experiment  are  shown  in 
Table  7.  This  variation  in  material  input  decision  shows  an  improve¬ 
ment  in  copper  production  rate  of  10.4  percent.  This  is  primarily  a 
result  of  the  smaller  number  of  conflicts  for  crane  service  in  run  2. 

Runs  3  and  4  illustrate  the  difference  in  performance  for  two 
different  initial  starting  conditions  as  given  in  Table  8.  No  signifi¬ 
cant  difference  in  performance  was  detected  for  these  particular  ini¬ 
tial  conditions,  however,  the  nunfcer  of  conflicts  for  crane  service  on 
each  converter  varied  considerably. 

Run  6  illustrates  the  performance  assuming  infinite  crane  service 
flexibility  (i.e.  the  crane  is  always  available  when  a  converter  requires 
service).  The  initial  conditions  are  as  in  run  3.  Table  9  shows  a 
comparison  between  run  3  and  run  6. 

A  production  rate  improvement  of  13.1  percent  is  obtained  by  assum¬ 
ing  infinite  crane  flexibility.  An  improvement  this  large  may  not  be 
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TABLE  6 

MATERIAL  INPUT  DECISIONS 


1 
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feasible  in  practice,  however,  because  a  no-conflict  crane  schedule 
require  increasing  blow  durations  in  some  cases  to  avoid  con¬ 
flicts. 


Optimization  Opportunities 

Summarizing  the  above  results,  we  observe  that  at  least  four  pos¬ 
sible  optimization  problems  may  be  identified.  These  include: 

1.  Optimum  crane  routing  so  as  to  minimize  the  time  a 
converter  must  wait  for  crane  service. 

2.  Optimum  scheduling  of  blow  durations  on  each  converter 
so  as  to  maximize  the  copper  production  rate  on  that 
converter. 

3.  Optimum  scheduling  of  activities  so  as  to  minimize  the 
number  of  conflicts  for  crane  service. 

4.  Control  of  the  converter  furnaces  in  accordance  with  blow 
duration  specified  by  the  crane  schedule. 

In  the  next  chapter,  these  problems  are  clearly  defined  and  sever¬ 
al  optimization  techniques  which  may  be  applied  are  presented  together 
with  some  simple  illustrative  examples. 


CHAPTER  V 


OPTIMIZATION 

If  it  were  not  for  the  mixed  continuous  control  variables  and 
discrete  decisions  with  unspecified  times  of  occurrence,  optimiza¬ 
tion  of  the  aisle  system  would  be  relatively  straight  forward.  As  is, 
the  problem  requires  modification  of  the  model  or  modified  applica¬ 
tion  of  standard  optimization  techniques.  In  the  following  sections 
we  focus  on  two  general  approaches  to  aisle  optimization.  First, 
for  direct  total  system  optimization,  continuous  variables  are 
quantized  in  time  and  combined  with  discrete  decision  variables  to 
form  one  large  discrete  time  model  which  is  theoretically  amenable 
to  solution;  for  the  second  approach,  the  primary  optimization  prob¬ 
lem  is  partitioned  into  a  number  of  smaller,  less  complex  optimization 
subtasks  and  these  are  related  to  total  system  optimization. 

Because  the  problem  for  the  direct  approach  is  so  large  and  comp¬ 
lex  that  it  is  unwieldy,  it  is  not  solved.  Instead  the  partitioned 
method  is  developed  to  the  point  where  a  solution  can  be  achieved. 

Direct  Optimization  Approach 

A  number  of  optimization  techniques  apply  to  large  discrete  time 
(or  discrete  state)  models.  The  aisle  model  given  in  Chapter  III  can 
be  converted  to  a  single  combined  model  with  this  form.  Let  A  be  a 
time  increment  such  that  time  duration  of  any  blow  or  any  crane  move¬ 
ment  can  be  represented  as  an  integral  multiple  of  A.  A  state  trans- 


ition  equation  for  the  ith  furnace  takes  the  general  form: 
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(5.1) 


where  j  =  0,  1,...  corresponds  to  t  =  0,  A,  2A,...;  x1  and  u1  are  the 

3 

ith  furnace  state  and  control  at  time  j ,  and  p ,  is  the  crane  state  at 
time  j  (which  defines  the  crane  service  at  time  j) ;  and  u  is  assumed 

j 


to  be  constant  for  jA£t_<(  j+1)  A . 


The  functions  gL  are  obtained  by  integrating  Equation  3.7,  and 
adjusting  the  state  to  account  for  possible  crane  service.  Of  course, 
Ej  only  affects  2£j+1  if  the  crane  services  furnace  i  at  time  j. 

Note  that  by  including  p  as  an  argument  in  g^,  Equation  5.1  models  the 
furnace  for  the  complete  cycle  including  the  times  of  the  crane  visit. 
Similarly ,  the  crane  transition  equation  can  be  written  as 


—j+1.  =  Jl(£j'  V  (5.2) 

where  p_.  and  d  are  the  crane  state  and  transition  decision  as  defined 
in  Chapter  III,  except  for  the  altered  choice  of  discrete  times.  Since 
the  discrete  times  are  not  restricted  to  the  times  of  the  crane  visits, 
(p^) j  assumes  the  discrete  states  0,  1,  2,  3  described  in  Chapter  III, 
and  an  additional  value  corresponding  to  "crane  in  transit."  Of  course 
the  crane  can  change  state  only  at  times  when  it  has  completed  a  pre¬ 
viously  initiated  move. 

Combining  Equations  5.1  and  5.2,  the  total  aisle  system  state 
transition  equation  takes  the  form: 


-j+1  ■  2(V  V 


(5.3) 
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and  £  is  given  as 
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The  system  constraints  and  the  performance  criterion  may  be  rec\dily 
expressed  in  terms  of  z_.  and  u .. .  Let  us  now  examine  the  magnitude  of 
the  optimization  task. 

The  model,  Equation  5. 3, still  has  mixed  control  variables  with 
continuous  and  discrete  values.  This  can  be  treated,  rut  for  simplicity 
assume  that  each  control  u_.  (oxygen  input  rate)  is  quantized  into  y 
values.  The  control  d.  has  five  pos£;ible  values;  but  these  are  not 
permissible  at  each  time  j;  hence,  assume  that  there  are  effectively  v 
crane  decisions  possible  at  each  time.  Assume  that  the  total  production 
cycle  lasts  N  time  intervals.  If  an  initial  state  were  specified,  there 
would  be  ( 3y  +  v)  possible  controls  a  number  which  easily  becomes 
unwieldly.  To  compound  matters,  the  initial  state  is  not  known.  Speci¬ 
fication  of  the  initial  state  is  equivalent  to  specification  of  the 
staggering  of  charge  cycles,  which  need  not  be  known. 

In  spite  of  the  enormity  of  the  problem,  direct  optimization  may 
be  feasible.  We  comment  on  three  techniques. ^ 

Dynamic  programming.-  This  is  a  highly  efficier *  technique  for 
solving  many  N-stage  decision  processes  [2].  The  technique  is  illustra¬ 
ted  later  using  a  much  smaller  process  model.  The  computational  sav¬ 
ing  afforded  by  dynamic  programming  is  that  the  number  of  controls  to 
be  evaluated  is  proportional  to  <3y  +  v)N  rather  than  (3y  +  v)N,  again 
assuming  an  initial  state  is  specified.  Larson  [6]  illustrates  a  fur¬ 
ther  computational  savings  over  direct  application  of  dynamic  programm¬ 
ing,  for  the  case  in  which  a  nominal  starting  state  can  be  specified. 


+_ 

Because  of  the  size  of  the  total  aisle  model,  no  examples  are 
presented  in  this  section,  but  optimization  techniques  are  illustrated 
m  the  partitioning  approach. 
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than  determining  the  optimal  control  for  each  possible  starting 
state/  only  initial  states  "near"  the  nominal  starting  state  are  examined. 

Gradient  search .  Gradient  search  in  the  space  of  control  sequen- 
ces  {v^,...,^}  is  normally  a  very  simple  technique  for  control  calcula¬ 
tion  [7,  10],  but  it  requires  that  the  system  performance  J  vary  con¬ 
tinuously  with  each  control  variable  (v  ) .  and  that  each  control  vector 

*  D 

Vj  be  independent  of  the  value  of  the  control  at  any  other  time. 

Unfortunately,  for  the  formulation  of  the  model  in  this  section,  the  con¬ 
trol  does  not  meet  these  requirements.  Thus,  the  gradient  technique  is 
not  directly  applicable  for  determination  of  the  optimal  aisle  control. 

On  the  other  hand,  the  gradient  technique  may  be  applicable  for  searching 
on  values  of  the  initial  state,  for  use  in  conjunction  with  some  other 

technique  for  determining  an  optimal  control  with  an  initial  state  speci¬ 
fied. 

Random  search. —  It  has  been  shown,  for  example  [4,14],  that  as  the 
number  of  decision  variables  in  a  problem  becomes  large,  a  random  search 
may  become  more  efficient  than  a  grad:  ent  search.  Intuitively,  the 
reason  for  this  is  the  gradient  must  be  computed  for  each  decision  vari¬ 
able  whereas  in  the  same  time,  the  random  search  is  likely  to  have  found 
an  improved  value  of  the  control.  Further,  control  constraints  do  not 
seriously  hamper  a  random  search  because  non-allowed  controls  are  merely 
ignored.  Thus,  random  search  may  be  a  candidate  for  direct  a'sle  optimi¬ 
zation.  Lest  the  random  approach  seem  unjustified,  we  note  that  it  is 
the  basis  for  EVOP  [3]  ,  an  "evolutionary  optimization"  scheme  which  has 
been  applied  to  control  of  chemical  processes. 
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Optimization  by  Partitioning 

Some  of  the  difficulties  encountered  in  the  direct  optimization 
approach  may  be  solved  by  use  of  the  principle  of  invariant  embedding 
[12] ,  according  to  which  a  very  difficult  or  unsolvable  problem  is 
embedded  into  a  class  of  simpler,  solvable  problems.  For  aisle  opti¬ 
mization,  we  consider  partitioning  the  problem  into  a  number  of  more 
simple  subproblems  which  relate  to  the  overall  objective.  The  struc¬ 
ture  of  the  aisle  problem  suggests  a  possible  hierarchy  of  control  sub¬ 
problems  as  listed  below: 

1.  Given  material  requirements  of  the  converters,  find  an 
optimum  crane  routing  which  will  minimize  the  service 
time  for  each  converter  service. 

2.  Given  material  input  decisions  for  a  converter,  find  the 
minimum  converter  blow  durations  to  minimize  the  total 
charge  cycle  time. 

3.  Given  fixed  crane  service  times  and  minimum  blow  dura¬ 
tions  for  each  converter,  find  an  optimum  converter 
schedule  which  will  minimize  the  charge  cycle  time  for 
each  converter. 

4.  Given  a  schedule  of  converter  services,  find  optimum 
continuous  inputs  (e.g.,  oxygen  input  rate)  to  each 
converter  furnace  which  will  achieve  that  schedule  and 
satisfy  the  physical  constraints  of  the  converting  pro¬ 
cess. 

A  systematic  structuring  of  these  four  problems  is  shown  in  Fig¬ 
ure  13.  Of  course,  this  structure  is  not  unique,  but  it  does  have  the 
following  advantages : 

1.  The  complexity  of  each  subproblem  is  less  than  that  of 
the  total  problem. 

2.  The  control  variables  in  each  subproblem  are  of  one  type 
(i.e.,  discrete  or  continuous). 

3.  The  subproblems  can  be  implemented  separately,  thus  allow¬ 
ing  opportunity  to  evaluate  the  results  of  one  problem 
before  proceeding  to  the  next. 


Assumed 
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The  primary  disadvantage  to  this  approach  comes  with  the  problem  of 
composing  a  total  system  solution  from  the  solutions  to  the  subproblems. 
In  the  fallowing  subsections  a  method  for  solving  each  of  the  above 

subproblems  is  presented.  Then,  in  Chapter  VI,  numerical  examples  are 
given. 

Crane  routing.—  Let  a  converter  aisle  material  handling  process 
be  characterized  by  the  network*"  shown  in  Figure  14.  The  nodes  in  the 
network  represent  crane  locations  and  the  numbers  on  the  branches  re¬ 
present  the  ladle  state  during  transit  and  the  transit  times  as  given 
in  Tables  10  and  11. 

The  crane  routing  problem  is :  given  requirements  for  converter 
services,  determine  the  sequence  of  crane  moves  to  meet  the  requirements 
in  minimum  time  for  each  converter.  Clearly,  if  only  one  converter 
needs  service,  the  crane  routing  is  trivial  and  the  service  time  is  im¬ 
mediately  determined.  If  several  converters  need  service  simultaneously, 
it  may  be  more  efficient  to  service  them  jointly  rather  than  sequentu- 
ally.  This  problem  arises  if  either  (i)  several  converters  have  been 
scheduled  for  simultaneous  service,  or  (ii)  in  actual  operation  of  an 

aisle,  random  fluctuations  cause  several  converters  to  need  services  at 
the  same  time. 

The  performance  criterion  of  the  total  system  was  given  by  Equation 
2.1  as 

j  -  I 

i=l  T. 

l 

+_.  .  .  . 

Figure  14  is  not  a  state- transition  diagram  because  the  states 
are  partially  associated  with  the  branches . 
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Figure  14  CRANE  ROUTE  NETWORK 
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where  the  maximum  of  J  is  sought  over  all  decisions  and  controls. 

Holding  all  decisions  and  controls  fixed  except  the  crane  routing 
decisions  d  ,  k  =  0,...,  K  required  to  meet  a  specified  service  re- 

X 

S 

quirement,  the  charge  cycle  time  may  be  decomposed  as  T\  =  T?  +  T\ 

g 

where  T|  is  fixed  and  is  the  time  required  for  the  service  to  be 

S  ; 

completed.  Because  is  small  compared  with  Tf,  the  performance  object¬ 
ive  can  be  approximated  as 


max  J  2  maximum  J  Cu. (T* 

L  i  i 

i=l 


<> 


do . dK 


which  is  equivalent  to 


max  J  ~  minimum 
d  ,  . . .  ,d 


3 

l 

i-1 


Cu.  T.  + 

l  l 


constant 


K 


This  determines  as  a  crane  routing  criterion  to  be  minimized 


JK  = 


K 

I 

k=0 


J,  Cui  ‘k  6ik 

i=l 


(5.4) 


where  t,  is  the  crane  transition  time  for  the  kth  crane  move  and  6 .  =1 

k  ik 

if  the  ith  converter  requires  service  during  the  kth  crane  move  and  6.  = 

ik 

0  otherwise. 

Formalizing  the  crane  routing  problem  now:  given  an  initial  crane 

state  and  initial  service  requirements  cy_nQ,  SLoutp'  determine  crane 

routing  decisions  d  ,...,d  such  that  q-  =  0,  qout  =  0,  p  =0  and 

J  is  minimum. 

K 

The  method  of  dynamic  programming  offers  a  systematic  approach  to 
problems  of  this  type.  The  essential  feature  of  this  method  is  given 
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by  Bellman's  "Principle  of  Optimality"  [1,2].  A  mathematical  statement 
of  the  principle  of  optimality  is  the  recursion  formula 


min  t,  +  J°  , 
,  k  k-1 

dk 


(5.5) 


where  J°  is  the  optimal  value  of  expressed  as  a  function  of  the  com¬ 
bined  state  (qin^/  Sout^'  •  in  practice,  the  problem  is  solved  using 

a  technique  of  successive  approximation.  It  is  assumed  that  the  optimal 
o 

Jk-1  1S  tabulated  for  each  combined  state  under  consideration  at  the 

(k-1) th  step.  For  each  state  at  the  kth  step,  J°  is  evaluated  by 

searching  over  all  tabulated  states  at  the  (k-l)th  step  and  all  values 

of  the  decision  dk.  The  procedure  eliminates  all  routes  to  each  state 

that  have  Jk  greater  than  minimum  for  that  state.  Incrementing  k,  the 

optimal  route  will  be  determined  if  the  final  state  is  achieved  for  all 

possible  routes  or  if  all  other  partial  routes  have  J°  greater  than  J° 

K  K 

for  a  route  which  achieves  the  final  state. 

To  illustrate  this  technique,  let  initial  vectors  be  =  [2  2  0]fc 

SLout0  “  [0  2  0]  ,  p^  =  [0  0]t.  Since  converter  three  is  not  involved, 
the  third  variable  will  be  dropped  from  the  queue  vectors.  The  dynamic 
programming  tabulation  of  J°  as  a  function  of  the  combined  state  is  shown 
in  Table  12,  where  the  eliminated  rows  were  nonoptimal  routes.  For  ease 
in  reading  the  table,  the  states  for  each  k  value  are  labeled  with  lower 
case  letters.  The  table  is  generated  up  to  k  =  10,  at  which  stage 
the  optimal  route  is  determined.  Tracing  backwards  through  the  table, 
the  optimal  sequence  of  states  is  marked  with  asterisks  . 
sequence  is 


The  optimal 


TABLE  12 


DYNAMIC  PROGRAMMING  CRANE  ROUTING  EXAMPLE 


p 

■ 

qink 

^out, 

k 

>1 

-1 

[0 

0] 

[2 

2] 

[0 

2] 

0 

[1 

1] 

[1 

2] 

[0 

2] 

10 

a 

* 

[2 

0] 

[2 

2] 

[0 

1] 

8 

a 

[1 

2] 

[0 

2] 

16 

a 

[1 

2] 

[0 

1] 

16 

a 

* 

[2 

2] 

[0 

1] 

14 

b 

[1 

1] 

[0 

2] 

[0 

2] 

26 

a 

[2 

0] 

[1 

2] 

[0 

11 

24 

a 

worse 

than  2b 

[0 

1] 

[1 

2] 

[0 

11 

22 

■ 

* 

[1 

1] 

[1 

2] 

[0 

1] 

24 

[2 

0] 

[2 

2] 

[0 

0] 

22 

[2 

0] 

[0 

2] 

[0 

1] 

29 

a 

[1 

1] 

[0 

2] 

[0 

1] 

32 

* 

[2 

0] 

U 

2] 

to 

0] 

30 

[0 

0] 

[1 

2] 

[0 

1] 

30 

B 

worse 

than  3c 

[2 

0] 

[1 

2] 

to 

0] 

30 

B 

same 

as  4c 

[0 

1] 

[2 

2] 

to 

0] 

28 

e 

[0 

11 

[0 

2] 

[0 

1] 

32 

a 

[2 

0] 

[0 

2] 

[0 

0] 

35 

b 

* 

[0 

1] 

[1 

2] 

[0 

0] 

36 

c 

[1 

1] 

[1 

2] 

to 

0] 

33 

f 

[2 

1] 

[2 

1] 

[0 

0] 

36 

f 

65 


TABLE  12  continued 


While  the  dynamic  programming  example  may  seem  complex ,  the  algo¬ 
rithm  can  be  implemented  simply  and  efficiently.  The  computational 
savings  over  direct  enumeration  of  all  possible  routings  is  evidenced 
by  the  early  termination  of  nine  partial  routes  in  Table  12.  By  com¬ 
parison,  only  one  route  which  reached  the  desired  state  was  evaluated. 

In  spite  of  these  arguments,  the  illustration,  example  was  contrived 
and  it  could  reasonably  have  been  solved  by  enumeration  of  routes.  On 
the  other  hand,  for  a  more  complex  aisle  with  additional  cranes,  con¬ 
verters,  and  reverberatory  and  holding  furnaces,  routing  problems  can 
become  less  trivial.  Also,  dynamic  programming  can  be  extended  to 
total  system  optimization  as  was  discussed  previously. 

Minimum  Blow  Durations 

Given  the  material  input  decision  for  each  blow  on  the  jth  convert¬ 
er,  the  problem  is  to  find  the  oxygen  input  rate,  u,  and  the  blow  dura- 

tions ,  TB (n)  which  maximize  _ j_  .  since  Cu  .  is  fixed  by  the  material 

T.  3 

3 

input  decisions,  then  maximizing  the  copper  production  rate  is  equiva¬ 
lent  to  minimizing 

n 

T  =  l  TB. (n)  (5.6) 

J  i=l 

where  TB (n)  is  the  duration  of  the  nth  biow.  If  u  is  to  be  held  constant 
at  the  maximum  rate  during  the  finish  blow,  then  the  duration  of  the 
finish  phase  blow  is  determined  by  the  relationship 
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TB,  (finish)  =  -  ln  X2(tf)  (5.7) 

1  u  — , 

X2(o) 

wner.  Jie  ratio  /y.^(o)  is  the  concentration  of  Cu^S  at  the  end 

of  the  finish  blow.  The  minimum  blow  duration  of  the  finish  blow  then 
is  specified  by  the  maximum  u  available.  Thus,  only  the  duration  of 
slag  blows  is  considered  in  the  optimization  problem. 

This  problem,  like  the  routing  problem,  yields  to  solution  by  us¬ 
ing  the  dynamic  programming  technique.  Expressing  the  state  equations 
for  the  slag  phase  in  terms  of  the  two  independent  state  variables  PeS 
and  temperature,  we  obtain 


where 


xi =  anxi +  aisx5 +  v  + 


X5  +  a51Xl  +  a55X5  +  V  +  C5 


(5.8) 


all 

i — 1 
«— 1 

li 

+  k3S13 

+  + 
4  14 

k6a!6 

a51 

a51 

+  k3a53 

+  k.a_,  + 
4  54 

k6a56 

C1  = 

C1  + 

al2X2 

C5  = 

C5  + 

a52X2 

where  unprimed  constants  and  variables  are  defined  in  Chapter  III. 
The  state  variables  are  quantized  into  q  levels  for  X  and  p  levels 
for  X2  where  and  H ^  vary  over  the  range 

0  <_  X  <_  max  tons 

(5.9) 

1900  <  X5  <  2250 °F 
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The  blow  time  decision  variable  TB  may  be  quantized  into  y  levels  in 
the  interval  given  by 


TB  min  £  TB  £  TB  max 

The  oxygen  input  rate  u  may  be  a  constant  during  the  slag  blow  or  a 

continuous  function  of  time.  If  u  is  a  continuous  function  the  search 

space  becomes  infinite  dimensional  and  thus  adds  a  higher  level  of 

complexity  to  the  problem.  For  simplicity,  assume  u  to  be  constant  dur- 

mg  each  slag  blow,  thus  it  may  be  quantized  into  m  levels  in  the  inter¬ 
val 


min  —  —  max 

Since  there  are  n  slag  blows,  the  problem  is  formulated  as  an  n  stage 

decision  process  as  illustrated  in  Figure  15.  The  state  transition 
rule  given  by 


-(V  =  1  X(n-l)  + 


rtfn 


^  b  u  d  t 


(5.10) 


where  X(tn)  is  the  state  vector  at  the  end  of  the  nth  blowing  period, 

x(n)  is  the  state  vector  at  the  end  of  crane  service  following  the  nth 

blow,  and  l  is  the  fundamental  solution  matrix  to  the  homogeneous  form 

of  Equation  5.8.  The  state  vector  at  the  end  of  crane  service  following 
the  nth  blow  is  given  by 


X(n)  X  (t^n)  +  Ax(e,m) 


(5.11) 


where  Ax  is  the  change  in  state  due  to  crane  service  in  Equation  3.3. 
At  each  stage,  all  possible  combinations  of  TB  and  u  are  considered 
for  each  of  the  states  attained  from  two  separate  states  at  stage  (n-1) 
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as  illustrated  m  Figure  16.  By  the  principle  of  optimality,  we 
choose  the  decisions  (TB,  u)  which  correspond  to  the  minimum  of 

Si  “  '"in  [Si-1  +  TB<“>]  (5.12) 

thus  trajectory  AC  is  the  optimum  path  to  state  C  at  stage  n.  Clearly 

each  of  the  states  at  stage  n-1  must  be  considered  for  all  decisions 
(TB,  u) . 

Converter  schedule.-  The  third  subproblem  is  that  of  several  con¬ 
verters  systematically  receiving  all  services  via  a  limited  facility, 

i.e.,  a  single  crane.  This,  leads  to  a  scheduling  problem  wherein  the 
objective  is  to  systematically  arrange  in  time  the  services  to  each  con¬ 
verter  so  that  the  total  production  rate  of  all  converters  is  maximized. 
Solution  of  the  scheduling  problem  results  in  a  schedule  for  the  crane 
operator,  which  specifies  the  times  crane  service  is  required  for  each 
converter,  and  a  set  of  schedules  for  the  converter  operators,  which  speci 
fy  the  blow  durations.  To  illustrate  these  ideas,  we  present  a  simple 
example  based  upon  the  following  assumptions: 

1*  1116  nuniber  blows  per  charge  cycle  is  given.  This 

determines  a  fixed  number  of  services,  one  between 
each  pair  of  consecutive  blows. 

2.  The  number  of  ladles  of  output  and  input  to  complete 
each  service  is  given. 

3.  For  simplicity,  it  is  assumed  that  once  a  service  is 
initiated,  it  will  be  completed,  i.e.,  all  output 
and  input  requirements  will  be  met,  before  service  to 
another  converter  is  begun,  m  other  words,  we  require 
a. schedule  with  no  conflicts  in  demands  for  crane  ser¬ 
vice.  As  previously  noted  service  times  are  known. 

assumption  3  may  be  unjustified.  The  crane  routing  example  showed  that 

it  can  be  more  efficient  to  service  two  converters  simultaneously  ra- 
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ther  than  causing  one  converter  to  wait  until  the  other  completes  ser¬ 
vice.  However,  the  routing  example  assumed  that  all  other  decisions 
were  fixed,  whereas  purposefully  scheduling  two  converters  for  simul¬ 
taneous  service  effects  all  other  times  for  service.  This  requires  fur¬ 
ther  investigation. 

With  the  above  assumptions,  generation  of  a  schedule  is  a  small 
enough  problem  that  it  may  be  solved  by  direct  enumeration  of  all  possi¬ 
ble  schedules  which  meet  all  the  system  constraints.  A  simple  example 
will  illustrate  the  scheduling  ideas.  Typical  service  times  and  mini¬ 
mum  blow  durations  for  a  three  converter,  chree-blow-per- charge  cycle  de¬ 
sign  are  shown  in  Table  13.  Service  times  are  denoted  as  TS(i,  j)  and 
minimum  blow  durations  by  TB(i ,  j) ,  where  i  denotes  the  ith  converter 
and  j  denotes  the  jth  blow.  Define  TO(i,  j)  as  the  scheduled  time  to 
service  the  ith  converter  prior  to  the  jth  blow.  The  system  constraints 
require  that 

TO (i '  3)  >  TO(i,  j  -  1)  +  TS(i,  j  -  1)  +  TB(i,  j  _  i)  (5.13) 

and 

TO(i ,  j)  i-  TO(n,  m)  +  a  TS(n,  m) 

for  n^i,  m  =  1,  2,  3,  and  0  £  a  £  1 .  The  objective  of  scheduling  is 
to  maximize  J,  Equation  2.1,  and  hence,  with  the  assumed  fixed  material 
input  and  output,  to  minimize  the  charge  cycle  times  T. ,  given  as 

3 

Ti  =  I  lTB(i,  j)  +  TS (i ,  j)]  (5.14) 

3=1 

where  TB(i,  j)  is  the  scheduled  blow  time,  which  is  determined  as 


TB(i,  j)  =  TO(i ,  j  +  1)  -  TO (i ,  j)  -  TS(i,  j) 


(5.15) 


TABLE  13 

SERVICE  TIMES  AND  MINIMUM  BLOW  DURATIONS 


Converter 

1 


Service 

Time 

TS(l,j) 

Minimum 

Blow 

Duration 
TB  (l,j) 

TS(2,j) 

TB  (2 , j ) 

TS(3,j) 

TB  (3 

40 

80 

37 

80 

39 

80 

30 

50 

25 

50 

29 

50 

28 

90 

22 

90 

25 

90 

Blow  j  2 


TABLE  14 

CRANE  SERVICE  SCHEDULE,  TO(i,  j) 


Converter 
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For  the  sample  data  in  Table  13,  a  schedule  was  determined  as  shown  in 
Table  14,  where  we  assume  T0(1,  1)  =0. 

In  order  to  achieve  the  blow  durations  determined  by  scheduling, 
we  turn  attention  to  control  of  a  converter  furnace. 

Converter  furnace  control. —  We  are  motivated  to  study  converter 
furnace  optimization  by  the  desire  to  control  the  blow  durations  in 
accordance  with  a  specified  schedule.  However,  there  are  several  other 
variables  that  might  be  of  interest  in  addition  to  blow  duration.  These 
include  temperature,  sulfur  dioxide  (SO  )  production  rate  and  magnetite 
formation.  Temperature  is  important  because  high  bath  temperatures  dam¬ 
age  the  refractory  lining  of  the  converter  and  at  low  temperatures,  mag¬ 
netite  is  formed,  which  is  generally  undesirable.  S02  production  rate 
is  important  since  efficient  operation  of  the  facilities  which  process 
S02  often  depend  on  a  constant  output  from  the  converter  aisle.  We 
new  survey  a  number  of  optirrud  control  solution  techniques. 

Perhaps  the  simplest  control  is  to  assume  that  the  oxygen  input 
rate  is  constant  over  the  entire  blowing  period.  The  behavior  of  the 
converter  furnace  is  determined  for  each  constant  value  of  oxygen  input 
by  solving  the  model  equations,  for  example,  Equations  3.9  and  3.13. 

For  a  linear  model  state  equation  of  the  form 

x(t)  =  ftx:(t)  +  B  u  +  C  (5.16) 

where  u  is  the  scalar  control  variable  (oxygen  input  rate) ,  the  solu¬ 
tion  takes  the  form 

t 

x(t)  =  _^(t,  tQ)  x(t  )  +  f  $_(t,  T)  [Bu  +  C]  dt 

•'t  ~ 

o 


(5.17) 
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where  £  is  the  fundamental  solution  matrix  corresponding  to  the  homo¬ 
geneous  system  [12] .  Since  we  have  assumed  u  to  be  constant,  it  can 
be  removed  from  under  the  integral  sign,  thus  yielding  an  equation  of 
the  form 


x(t)  =  y^t)  +  y.2(t)  u 


(5.18) 


where  y  (t)  and  y2(t)  are  known,  Substituting  various  values  of  u  into 
the  above  equation,  we  obtain  state  trajectories  which  are  to  be  evalu¬ 
ated.  By  choosing  from  these  the  one  that  represents  the  most  desirable 
behavior,  we  have  in  a  very  simple  way  solved  the  problem. 

If  a  constant  control  is  not  adequate,  or  it  does  not  satisfy  the 
constraints  (Equation  3.13),  then  we  define  a  performance  criterion  J 
and  determine,  in  general,  a  time  varying  control  function  which  will  . 

/v 

influence  the  converter  to  minimize  J.  To  illustrate,  we  pose  the  fol¬ 
lowing  problem:  using  the  slag  blow  model  developed  in  Chapter  III  and 


assuming  SC>2  production  rate  to  be  linearly  proportional  to  FeS  oxida¬ 
tion  rate,  we  express  a  performance  criterion  as 


J  =  min 
u 


* 


<Vi 


R)  2  dt 


(5.19) 


where 


t^  =  scheduled  blow  duration 
K^x.,  =  S02  discharge  rate 
R  =  desired  SC2  discharge  rate 
u  =  oxygen  input  rate 


We  offer  as  a  justification  for  this  performance  criterion,  the  physi- 
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cal  argument  that  it  is  desirable  to  maintain  a  constant  so.,  discharge 

during  each  blow  duration.  The  rate  of  oxidation  of  FeS  was  expressed 
by  the  model  Equation  3.9  as 


*1  -a1*+b1u+  c2  (5. 20) 

where  a^  -  [a^  •••a^].  Substituting  Equation  5.20  into  Equation  5.19 
and  combining  constants  yields  J  of  the  form 


J  =  min 


(5.21) 


Note  that  J  is  a  quadratic  function  of  u  and  thus  u  will  be  directly 
determined  by  .minimizing  J. 

The  system  constraints  are  rewritten  as 

^  ^  xi,  /  k  =  1,  5  /  c. 

K  f  —  k  max  '  (5.2x 

0  1  *k ( tf ) ,  k  =  2 ,  3 ,  4 ,  6 


0  <  u  <  u 

-  ~  max  (5.23) 

We  constrain  the  terminal  value  of  the  state  variables  via  Equation  5.22 
rather  than  constraining  the  state  trajectories,  because  of  the  mono¬ 
tonic  nature  of  the  state  variables  and  because  the  terminal  constraint 
problem  is  generally  easier  to  solve. 

With  the  problem  now  expressed  mathematically,  we  point  out  two 
dynamic  optimization  techniques:  (  )  solution  via  a  set  of  necessary 
conditions  determined  by  the  maximum  principle  [7,  10],  and  (ii)  solu¬ 
tion  via  a  gradient  search  in  the  space  of  control  functions  [7,  10] . 

Both  techniques  employ  a  function  II  called  the  Hamiltonian,  which  is 
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based  upon  calculus  of  variations. 

For  the  approach  via  the  maximum  principle,  the  constraint  Equa¬ 
tions  5.22,  5.23  are  absorbed  in  two  ways,  and  then  the  set  of  neces¬ 
sary  equations  are  written  which  specify  the  optimal  u°(t).  Equation 

5.22  is  absorbed  by  attaching  a  penalty  function  [7,  10]  to  Equation 
5.21. 

A 

J1  =  3  +  WlXi(tf)  +  wiW  (5.24) 

where  and  w5  are  penalty  weights  which  must  be  determined,  perhaps  on 

a  trial  basis,  so  that  Equation  5.22  will  be  satisfied.  Equation  5.23 

is  absorbed  by  defining  new  variables  y  and  y  ,  where  [y  (t) ] 2  =  u(t)  , 

■*■2  1  f 

and  using  the  new  constraint 

Y1  +  Y2  =  Umax  (5-25) 

The  Hamiltonian  for  this  problem  is 

H  =  (“i*  +  +  ^  +  A  (Ax  +  By ^  +  C)  +  T(y^  +  y2  -  u  )  (5.26) 

—  1  ~  12  max 

where  r  is  a  scalar  Lagrange  multiplier  and  A(t)  is  a  vector  LagrcJige 
multiplier.  The  necessary  equations  are 


3H 

dX 


x 


(5.27) 


3H 

3x 


-A 


subject  to  the  boundary  conditions 
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( t^.)  “  k  =  1 ,  5 ;  X  =  o ,  k  =  2  3  4  5 

wnere  ^  is  assumed  given  and  and  w5  must  be  determined.  Equations 
5.27  together  with  the  boundary  conditions,  constitute  a  two-point 
boundary  value  problem  [7,  JO]  which  yi,  -  die  optimal  control  u°(t) . 

In  general,  such  a  system  is  difficult  to  solve,  but  for  this  particular 
case,  Equation  5.27  reduces  to  a  linear  differential  system  which  may 
be  readily  solved. 

In  some  cases,  control  solution  via  the  above  approach  may  be  too 
difficult.  For  example,  if  -.he  performance  criteria  are  not  a  quadratic 
function  of  the  control  u,  the  control  may  be  singular  [8].  In  these 
cases,  the  control  may  be  determined  by  an  iterative  gradient  search, 
which  is  simpler  but  less  direct.  At  each  iteration  of  the  gradient 
search,  an  approximate  u(t)  is  known  and 

3H  _  ,o 
3x  - 


are  solved  subject  to  known  boundary  conditions.  Then,  ~  is  computed 
and  used  to  update  the  approximate  u(t) ,  and  the  next  iteration  is  be¬ 
gun.  Nononen  and  Paquurek  [8]  illustrate  a  gradient  search  technique 
using  flux  addition  rate  as  the  continuous  control  variable  and  tempera¬ 
ture  as  a  controlled  variable. 

Thus  far  we  have  assumed  the  blow  duration  to  be  fixed  by  a  sched¬ 
ule,  based  upon  an  estimate  of  the  minimum  blow  times.  If  the  scheduled 
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blow  times  were  not  realistic,  however,  it  may  not  be  possible  to 
find  a  satisfactory  solution,  or  any  solution  for  that  matter,  with 
the  blow  terminal  time  fixed.  This  is  motivation  to  consider  a 
variable  end  time  problem.  Two  basic  methods  of  dealing  with  this  prob¬ 
lem  are:  (i)  to  imbed  time  into  the  original  problem  as  a  new  variaole 
and  optimize  using  one  of  the  above  techniques,  or  (ii)  simply  to  search 
on  the  end-time  by  solving  the  problem  for  various  values  of  ending 
time.  Solution  of  the  variable  end-time  problem  may  result  in  blow 
durations  which  are  inconsistent  with  the  schedule.  To  resolve  this 
difficulty,  we  turn  to  the  final  task  of  coordinating  the  subproblems 
that  have  been  presented  so  as  to  achieve  a  total  system  optimization. 


CHAPTER  VI 


OPTIMIZATION  RESULTS 

In  the  previous  chapter,  several  optimization  problems  are  pre¬ 
sented  along  with  a  survey  of  techniques  for  solving  them.  We  now 
focus  attention  on  the  optimization  of  the  overall  converter  aisle 
performance.  To  achieve  the  goal  of  system  optimization,  the  problem 
is  partitioned  into  three  subproblems  as  shown  in  Figure  17.  The  crane 
routing  subproblem  is  not  included  in  this  study  since  it  applies  prim¬ 
arily  to  converter  aisles  having  a  large  number  of  converters  say01-  t,<- 

r  1  eignt  or 

nine,  while  this  is  an  example  of  a  three-converter  aisle.  To  eliminate 
this  problem,  we  schedule  to  avoid  crane  conflicts.  Assuming  material  input 
decisions,  the  minimum  blow  durations  and  corresponding  terminal  states 
are  determined.  Using  these  minimum  blow  durations  and  materia.1  input 
decisions,  an  optimum  no-conflict  schedule  is  calculated.  As  a  result 
of  scheduling,  the  blow  duration  may  be  longer  than  the  minimum,  thus 
the  converter  furnaces  are  regulated  in  accordance  with  these  durations 
by  manipulating  the  oxygen  input  rates.  The  important  advantage  of  the 
partitioned  approach  is  that  the  problem  has  been  reduced  to  a  search 
problem  on  a  small  set  of  decision  variables,  namely,  the  material 
input  decisions.  The  subproblems  may  be  formulated  as  standard  optimi¬ 
zation  problems  and  yield  to  solution  by  well- developed  techniques. 

For  simplicity  in  solving  each  of  the  subproblems  of  Figure  17,  the 
following  assumptions  are  made  held  constant  at  any  level  within  a  g  ven 


range . 
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1-  The  oxygen  input  rate  is  constant  during  each  blow. 

2.  Material  inputs  are  identical  for  all  converters. 

3.  Five  slag  blow  are  assumed. 

4.  The  system  is  not  subject  to  random  perterbations . 

Using  these  assumptions,  numerical  examples  of  the  three  subproblems 
of  Figure  17  are  presented.  Finally,  the  converter  aisle  is  optimized 
for  three  choices  of  material  inputs  and  the  results  are  compared,  but 
a  total  search  of  the  material  input  space  is  not  done. 


Minimum  Blow  Durations 

The  objective  is,  given  material  input  decisions,  to  find  the  oxygen 
input  rate,  u,  and  minimum  blow  durations  TB (n)  which  minimize  the 
objective  function  of  Equation  5.6. 

Since  it  is  necessary  to  numerically  integrate  Equation  5.11  for 
each  investigation,  the  computer  time  required  to  solve  the  problem 
(assuming  one  second  per  integration)  could  be  several  hours. 

To  establish  the  quantization  size  for  the  states  and  controls,  it  is 
necessary  to  try  several  combinations  of  each.  To  shorten  the  time 
required  to  obtain  a  solution,  the  state  model  is  approximated  by  the 
linear  relationship 


Xl(tn}  =  xX(n  "  +  -00177  u  TB (n) 

X5  (tjj)  =  X5  +  -0156  u  TB(n) 


(6.1) 


This  approximation  is  reasonable  for  this  problem  since  the  model  solu¬ 
tions  shown  in  Figure  11  are  linear  for  Xj_  and  x,..  An  analysis  of  the 

solutions  to  Equation  5.10  for  various  values  of  u  shows  x  and  x  to 

1  5 
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be  linear  in  u.  To  complete  the  state  transition  model  for  a  given 

state  Equation,  5.11  is  used.  The  amount  of  Cu  S  (i.e.  x  )  does  not 

2  2 

change  during  a  slag  blow  due  to  blowing  but  accumulates  with  each 
material  addition.  Thus  the  state  transition  equations  of  6.1  do 
not  reflect  changes  in  x2  but  Equation  5.11  does.  The  solution  begins 
at  the  initial  state  and  solves  forward  by  considering  all  combinations 
of  TB  and  u  for  all  states  x  at  each  state.  The  solution  has  the  con¬ 
straint  that  all  FeS  must  be  oxidized  at  the  end  of  the  fifth  blow. 

The  principle  of  optimality  given  by  Equation  5.12  specifies  the  decision 
rule  for  selecting  the  optimum  decisions  on  TB  and  u.  The  dynamic  pro¬ 
gramming  procedure  is  implemented  on  a  digital  computer  and  a  listing 
of  the  program  is  shown  in  Appendix  3.  The  results  of  a  dynamic  pro¬ 
gramming  solution  using  the  material  inputs  of  Table  6  Run  2,  are  given 
in  Table  15.  For  this  solution,  the  states  are  quantized  into  5  ton 

increments  for  0  _<  x  _<  30  tons  and  50  °F  increments  for  1900  °F  <  x  < 

x  —  2  — 

2250  F.  The  controls  are  quantized  into  10  minute  increments  for  30 
min  TB  jc  100  min  and  50  lb/minute  increments  for  100  <  u  <  250. 

Scheduling 

Given  the  material  inputs  of  Table  6,  Run  2; the  crane  service 
times  TS (i , j)  as  given  in  Table  16;  and  the  minimum  blow  durations  s 
given  in  Table  15;  find  a  schedule  of  converter  aisle  activities 
which  will  maximize  the  copper  production  rate  given  by  Equation  2.1. 

3  Cu. 

J  =  l 

i=l  T . 
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TABLE  15 

MINIMUM  BLOW  TIMES 
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Cu^  is  a  constant  determined  from  the  material  input  vector  by  the 
relationship 

6 

Cu  =  l  m  GW  (6.2) 

j=l  3 

where  G  =  matte  grade  (percent  Cu) 

W  =  weight  of  matte  per  ladle 
m  =  number  of  ladles  of  input  material 
Since  Cu^  is  a  constant,  the  maximization  of  Equation  2.1  is  achieved 
by  minimizing  the  charge  cycle  time  T\  which  is  given  by  equation 
5.14 

6  _ 

T.  =  l  [TB(i,j)  +  TS (i ,  j)  | 

The  solution  is  a  schedule  of  times  to  start  the  ith  blow  on  the  jth 
converter  and  a  Table  of  converter  blow  durations.  The  solution  is 
subject  to  the  constraints  of  Equation  5.18.  A  computer  program 
which  performs  an  exhaustive  search  on  all  possible  schedules  is 
shown  in  Appendix  4.  On  each  iteration  of  the  search,  the  value  of 
J  is  computed  and  compared  against  the  previous  value.  The  schedule 
associated  with  the  optimum  choice  is  retained  and  another  iteration 
is  begun.  The  starting  time  T0(1,1)  is  assumed  to  be  zero  in  every 
iteration.  The  results  of  this  procedure  using  the  material  inputs 
of  Run  2  from  Table  6,  the  minimum  blow  durations  from  Table  1^  and 
the  crane  service  times  for  Run  2,  Table  16  are  shewn  in  Table  17  and 


and  18. 
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Having  specified  the  blow  durations  as  given  in  Table  18,  we 
turn  to  designing  the  oxygen  input  rates  which  will  control  the 
oxidation  rates  of  the  converters  in  accordance  with  these  durations . 


Converter  Furnace  Control 

Regulation  of  the  converter  furnace  activity  to  achieve  a  fixed 
blow  duration  is  a  terminal  control  problem.  In  general,  the  object¬ 
ive  for  the  slag  blows  is  to  find  a  continuous  control  u°(t)  which 

transfers  the  system  from  an  initial  state  x(t  )  =  x  to  a  terminal 

—  o  — o 

state  given  by 


0  1  V  -  ^kmax  k  =  ifS  (6.3) 

in  a  time  interval  [tf  -  tQ]  =  TB(i,j)  for  j  -  1 . 5.  For  the  finish 

blow,  the  objective  is  to  transfer  the  initial  state  x  (t  )  =  X  to 

2  o  20 

the  terminal  state  x^ (t^)  *  0  in  the  interval  [t^  -  t  ]  =  TB(i,6) 

The  type  of  control  functions  u(t)  which  is  to  be  obtained  depends 

on  the  method  of  its  implementation.  If  the  control  adjustments  are  to 

be  made  manually  by  an  operator,  they  must  be  much  less  frequent  than 

the  case  where  an  automatic  controller  is  used.  Thus,  the  optimization 

method  for  finding  the  control  is  chosen  in  accordance  with  the  desired 

form  of  the  control.  For  manual  adjustment  u  must  be  a  discrete  form 

where  the  period  between  adjustments  is  sufficient  to  allow  operator 

response.  For  automatic  control,  u  may  be  a  continuous  function  of 

the  form  u  =  u(t)  for  all  t  in  [t  ,  t  ] . 

o  f 

In  Chapter  V,  two  techniques  for  determining  the  optimum  continu¬ 
ous  oxygen  input  rate  are  presented.  The  Maximum  Principle  of  Pcntry- 
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agin  [10]  specified  the  necessary  conditions  for  an  optimum  solution. 
The  gradient  search  technique  was  considered  as  an  alternative  method 
for  a  problem  where  the  control  is  singular  (i.e.,  the  control  cannot 
be  expressed  explicitly  as  a  function  of  the  state  variables  and 
Lagrange  multipliers) . 

For  a  discrete  solution,  the  simplest  control  is  to  assume  the 
oxygen  input  rate  is  constant  over  the  entire  blowing  period.  As  a 
practical  solution,  the  constant  oxygen  input  problem  is  considered. 

The  problem  is  to  find  the  constant  oxygen  input  rate  which  will 
transfer  the  stal  e  from  a  given  initial  value  to  a  given  terminal 
value  in  a  giver,  period  of  time  TB  (n)  ,  subject  to  the  constraints  of 
Equation  5.9.  The  oxygen  input  rate  is  implicitly  constrained  by  the 
solution  obtained  in  the  minimum  blow  time  problem.  An  iterative 
search  technique  is  used  to  find  u.  This  is  accomplished  by  solving 
Equation  5.10  for  various  values  of  u  and  testing  to  determine  if  the 
terminal  constraints  have  been  achieved.  A  computer  program  which 

performs  an  iterative  search  for  the  optimum,  u,  is  shown  in  Appendix 
5. 


The  control  for  the  finish  phase  may  be  obtained  directly  by 
solving  the  finish  phase  state  equation  and  inverting  thus 


u 


In  X2(tf} 


X2(o) 


(6.4) 


where  X2(tf)/X2(o)  is  the  terminal  concentration  of  Cu^S. 
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A  numerical  example  for  a  solution  to  this  problem  is  show,  in 
Table  19.  The  blow  durations  TB<„)  are  obtained  from  Table  la  and 
the  initial  and  terminal  states  are  obtained  from  Table  15. 

Comparative  Results 

Three  optimization  studies  were  concocted,  corresponding  to 

material  input  decision  of  Run  1  and  Run  2  in  Table  6  and  Run  3  of 
Table  20. 

Applying  the  optimization  technique  for  finding  minimum  blow  dur¬ 
ations  to  the  material  inputs  of  Run  1  yield  a  solution  which  does 
not  satisfy  the  constraints  of  the  minimum  blow  time  optimization 
problem,  namely  that  all  Res  must  be  oxidized  at  the  end  of  Blow  5. 
Therefore,  Run  1  is  eliminated  from  future  consideration  as  a  material 
input  choice.  In  practice,  when  the  situation  arises  where  Res  remains 
in  the  bath  at  maximum  temperature,  cold  copper  scrap  is  added  to 
cool  the  bath  and  blowing  is  resumed  until  all  ReS  has  been  oxidized. 

The  results  of  optimizing  for  the  material  inputs  „f  Runs  2  and 
3  are  compared  in  Table  21.  comparing  the  optimized  Run  2  in  Table 
21  With  the  non-optimized  R,„  2  in  Table  7  shows  an  improvement  of  29 
percent  in  favor  of  the  optimally  controlled  aisle.  Only  a  slight 

difference  is  noted  between  the  results  of  Runs  2  and  3  for  the  opti- 
mally  controlled  aisle. 

we  now  investigate  the  size  of  the  material  input  search  space. 

The  range  of  material  input  choicer  f nv 

put  cnoices  for  this  example  is  constrained  by 

the  physical  size  of  the  converters  and  the  assumption  that  only  an 

integer  number  of  ladles  may  be  transferred  during  crane  service. 


1 


92 


TABLE  20 

MATERIAL  INPUT  DECISIONS  FOR 
RUN  3  OF  THE  OPTIMIZATION  STUDIES 


Blow 


Number  of  Ladles 
_ of  Input _ 


4 

4 


93 


CM 

W 

-J 


t* 

o 

w 

H 

2 

H 

to 


to 
H 
P 
O  P* 

z  z 


H 

3 


►J 

.  <! 
W  M 
P<  Pi 

o  w 

H 

lg 

M  O 
H  W 
P<  XI 
O  H 


ft, 

o 

CO 

H 

g 

CO 

w 

&< 


X 

X 

H 


94 


Assuming  the  maximum  capacity  of  the  furnaces  to  be  M  tons,  the 

c 

objective  is  to  determine  the  number  of  ladles  of  material  and  the 
number  of  blows  which  minimize 

i 

M  -  l  MWG  >0  (6.7) 

j=l  11 

where  W  =  weight  per  ladle  of  matte 

=  weight  fraction  of  Cu2S  per  ton  of  matte 
i  =  number  of  slag  blows 
ML  =  number  of  ladles  on  the  jth  blow 

For  this  example  the  range  of  i  is  3  to  13  and  the  range  of  M  is  1 
to  9. 

Because  of  the  accumulation  of  Cu2S  on  successive  blows  and  the 
capacity  constraints,  the  search  space  is  approximately  200  points. 
With  this  size  of  space,  direct  enumeration  is  possible.  However,  a 
search  may  be  more  efficient. 
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CONCLUSIONS 

A  complex  industrial  optimization  problem  exemplified  by 
the  converter  aisle  in  a  copper  smelter  has  been  identified. 

Two  approaches  to  solving  this  optimization  problem  have  been 
presented.  The  direct  approach  offered  well  developed  methods 
for  solution,  but  the  problem  was  too  large  from  a  practical 
point  of  view.  Thus,  an  alternative  method  of  optimization  by 
partitioning  the  problem  into  managable  subproblems  was  developed. 
Then  the  solutions  to  the  subproblems  together  compose  a  solu¬ 
tion  to  the  total  problem.  The  results  of  applying  the  parti¬ 
tioned  optimization  technique  showed  that  an  improvement  of  at  least 
20%  in  the  copper  production  rate  was  achieved.  This  was  accomp¬ 
lished  by  minimizing  the  blowing  times,  scheduling  the  crane  ser¬ 
vice  activities, and  regulating  the  converter  furnaces  in  accord¬ 
ance  with  an  optimum  schedule.  Some  comments  are  in  order  re¬ 
garding  the  interpretation  of  the  findings  in  the  numerical  example. 

The  model  used  throughout  this  work  is  hypothetical.  Its 
purpose  is  to  characterize  the  dynamics  of  a  converter  aisle  with 
sufficient  realism  for  posing  a  practical  optimization  problem. 

The  important  contribution  of  this  research  is  the  development 
of  a  technique  for  solving  this  type  of  optimization  problem, 
which  otherwise  had  not  yet  been  solved.  To  view  in  perspec¬ 
tive  the  meaningful  results  of  the  partitioning  approach, 
two  important  sensitivity  considerations  should  be  recog- 
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nized.  First,  consider  the  sensitivity  of  the  model  parameters 
to  actual  process  operating  data.  The  form  of  the  model,  for 
example,  may  significantly  influence  the  sensitivity  of  the  model 
to  the  data  used  in  its  formulation.  One  may  also  ask, 

"How  are  the  results  obtained  from  the  model  changed,  by  varying 
the  model  parameters  (e.g.,  the  coefficients  in  Equation  3.9)?" 

This  question  has  not  been  investigated  in  this  work  but  clearly 
is  a  prerequisite  for  accepting  the  results  of  any  specific  study 
where  this  optimization  method  is  applied.  A  second  consideration 
is  the  sensitivity  of  the  optimization  results  to  the  model  of  the 
process  being  optimized.  Future  work  on  this  problem  should  include 
an  investigation  of  results  obtained  by  optimizing  with  several 
models. 

Some  considerations  which  may  influence  the  results  obtained 
from  the  optimization  approach  are: 

1.  The  numerical  integration  technique  used  to  solve  the 
model  equations . 

2.  The  step  size  of  the  numerical  integration. 

3.  The  quantization  levels  used  for  the  state  and  control 
variables  in  the  minimum  blow  duration  problem. 

Although  these  considerations  were  not  investigated  in  this  paper, 
it  should  be  recognized  that  they  may  have  an  important  influence 

on  the  results  obtained  from  an  actual  application  of  the  optimiza¬ 
tion  technique. 

Practical  application  of  this  technique  requires  initially  an  abun¬ 
dance  of  data  from  which  to  build  and  verify  a  model.  For  many  pro¬ 
cesses  such  as  the  converter  aisle,  the  type  of  data  needed  to 


accomplish  this  is  very  difficult  to  obtain,  principally  due  to  the 
lack  of  instrumentation  for  measuring  the  state  variables.  It  should 
be  recognized  that  the  physical  environment  is  very  harsh  and  special 
instrumentation  is  required  for  measuring  many  of  the  variables.  Cnee 
a  model  has  been  developed  and  verified,  application  of  the  optimiza¬ 
tion  techniques  presented  herein  will  provide  an  optimal  operating 
strategy.  However,  a  strategy  such  as  the  one  given  in  Table  22  does 
not  account  for  perturbations  in  the  process  such  as  changes  in  matte 
9^.ade ,  oxygen  input  rates,  or  equipment  failures.  To  compensate  for 
this,  we  seek  closed  loop  (i.e.  feedback)  control  which  continuously 
monitors  the  activities  of  the  process  and  adjusts  the  control  variables 
in  accordance  with  the  optimum  operating  strategy.  This  eventually 
requires  a  control  computer  which  performs  the  functions  of  monitoring 
the  system  variables,  filtering  measurements,  estimating  the  system 
state,  and  computing  the  optimal  control  values.  The  necessity  for 
feedback  control  may  be  a  factor  in  rejecting  the  direct  optimization 
approach  as  a  method  of  solving  this  problem.  In  the  case  where  the 
process  disturbances  are  such  that  it  is  necessary  to  compute  new 
optimal  control  values,  the  direct  approach  would  take  a  prohibitive 
length  of  time  to  compute  an  updated  operating  strategy  while  the  par- 
titioned  approach  produces  a  solution  in  a  much  shorter  time.  The 
direct  approach  does  have  an  advantage,  however,  in  that  a  true  optimal 
solution  is  guaranteed.  For  the  partitioned  approach  that  claim  cannot 
be  made  since  the  composition  of  the  subproblems  into  a  total  solution 
may  result  in  a  suboptimal  control  which  is  a  local  optimum. 
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Further  research  might  include  continued  study  of  the  direct 
approach  to  seek  methods  of  reducing  size  and  solution  time.  Further 
research  should  definitely  be  done  on  the  partitioned  approach  to 
consider  improved  models  and  the  significance  of  the  simplifying 
assumptions . 
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Appendix  1 

Computer  Program  of  the  Simulation  Model 
Subroutine  (EVNTS) 


SUBROUTINE  £VN T SI  I  X, NSFT ) 

0  IMENS  ION  N  SET ( 6 •  1  ) 

COMMON  10.  IM, INIT.JCVNT, JMNI T.MFA , MS  TO * , MX , MXC . NCLCT . NHIS T 

inoo.nikpt.n'jt.nprms.nrun.nr  uns  .nsta  t.out.scale.iseed  tnow 
ptbeg.tfin.mxx.nprnt.ncror.nep.vnoia)  •SCA‘-E.ISEED.TN0W. 

COMMON  ATR 13(4) ,ENO( 41  ,INN( 4), JCELSI5.2?)  .KRANKI4 I , J.  L R . M AX NO, 4  I 
IMFEI  4)  ,MLC(  4),  ML  El  4  )  .NCSLSI  51  ,N0<  4i  ,PAUAM(  ;«0 ,4  )  .  QT  I  ME  I  4  >  ,  s  SOMA  1 1  0 
25)«SUMA(  1  0 , 5  )  «  N  AME  I  6  )  «  NPRO  J  «  MCIN  «NOA  V  ,  NTR  * 

I  ^DELT^CU  ('3')"  ‘  ®  ^  ,  CQN0  ,  I  C  ,  I  O  «  NB  .TCU  (3)  .TSI3.6) 

I E=ATR I G { 4 ) 

I C  =  ATR  I B ( 3 ) 

V  I C=  IC 

CALL  F INO ( V  I C. 5 • 2, 2, KCOL .NSE T ) 

CALL  R  MOVE ( K  COL . 2, NSFT ) 

V  AL=  ATP  IB<  I  ) 

GO  Toll,?),  Ix 

IF (CPNP-I . ) 3,4,  4 
CALL  CRNSV(TSV) 

CALL  CCLCTI VAL .  IC.NStT ) 

ATRIBI  I)=TSV4TN0W 
I F (  IQ-MC 15.6,6 

ie=  is* 1 

GO  TO  7 
IB=  I 

TCUI  IC)=TCU(  IC  >  4CUI  1C) 

ATR IB( 4 )-  IH 
AT  R  IB( ?)=?, 

ATR  IB(  3 )= ic 
CRNE  = I . 


SAV= ATR IB( I )-TNOW 
CALL  FILEMI 1 ,NSET) 

ATR  10(1  )  =  SAV 

ATRIBI?)- IC 

CALL  F  ILFM<  ?,NSET ) 

RETURN 

A  CALL  FIND! TNU  W  «?,  I,  I,K  COL , N  SE  T ) 
CALc  RMOVE (KCOL ,  I  , NSFT ) 

« TM=ATP IB( I l-TNOW* .2 
IC4= IC43 

call  COLCTIWTM, ICW.NSET) 

CALL  TMSTI CRNfl, TNOW, I ,NSET) 

X1A  IT- ATR  I B (  1)4.? 

CALL  F  IL EM (  I , N  SE T ) 

ATR  IB(  I  >  =  WA  IT 
ATR  IB( 4  )=  13 
ATRIBI 3)= IC 
ATR  10(2)= I , 

CALL  F  ILEM (  I ,n SET ) 

ATR  IB  I  I )  =  VAL 4WTM 

ATR IB  I  2  >=  IC 

CALL  c  IL  EM  (2, NSFT) 

RETURN 
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Appendix  1  (Cont'd 


Subroutine  (EVNTS)  Cont'd. 


2  ATPI0(4»=I0 

ICP= IC+6 

CALL  COLCT ( V AL •  I  C°  *N  SET  J 
CALL  CONVT ( TMC I 
ATRI0C 1 >=TMC 
PAV=TNO**ATR IB<  1  ) 

CALL  F IL  EM ( 2  *N  SET ) 

ATP  I0« 1 JsPAV 
ATRI0( 2>= 1 . 

ATRI0(3)= IC 
ATR I0( 4 )=  |0 
CRN0=O . 

CALL  F  ILEM (  l.NSET) 

RETURN 

SNC 


Appendix  1  ( Cont' d.  ) 
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Subroutine  (CRfISV) 


♦ONE  WORC  INTEGERS 

SUBROUTINE  CRNSVITSV) 

01  MENS  IUN  PH  3) 

COMMON  in,  IM,  INIT.JEVNT, JMNI T , MF A , MS  TOP . MX , MXC , NCLC T  ,  NHI  ST  , 

1  NflO^  NC1RPT  »  NOT,  NPRM  S  ,  NR  UN  ,NR  UNS  ,N  STA  T  ,  OUT  (SCALE  ,  I  SEED  ,T  NOW  , 

2T0EG, TF  IN, MX  X , NPRN T , NCPOR ,NEP , VNOI  4  ) 

COMMON  ATRIUI4), ENO 14)  ,  I  NN< 4  I , JCt LS ( 5 . 22 > , KRANK ( 4 ) , JCLR , MAX NOI 4 )  , 
1MFEI 4  )  ,MLC( 4  ).MLEI 4)  .NCELSI  5) .NCK  V) .PAR AMI  2 0,4) , OT t ME! 4 ) , SSUMA (  1 0  , 
25 ) , SUM A(  1 0. 5 ) .NAME I  6  )  ,NPRO J ,MON ,NDA Y,NVR 

COMMON  XMATLI 1.6). X( 3,6) ,XN{ 14) , XD ( 6 ) .CRNB.IC . I B , NB , TCU (3 ) ,T5 I 3. 6 ) 
I. CELT,  CUI31.0XYI  t) 

TPLI  1  )  =  .7 
TPL( ?)= .5 
TPLI J)= . 7 
1 F I  IB-614, 5 , 4 
5  XL  CL  =  CU  I  10/15. *1. 

YLCL-O. 


GO  TO  7 

4  XL  CL- X (  IC, 4)/15,*l . 

VLCL=XL0L-1 . 

7  LCL- XL  CL 
XLCL  =  LCL 


T  SL -  XL  CL*  TPL I  1C  i 
1F<  10-5)1,2,  1 

C  COMPUTE  CONVERTER  STATE  AT  THE  END  OF  CPANE  SERVICE  FOR  SLAG’  HL  OW 
1  X  I  1C,  1  )  =  XMATLI  1  ,  10  )♦  XI  IC  ,  1  ) 


XI  IC, ?)= XMATLI 2,  1  3  )  ♦  X|  1C, 2) 
XI  IC,1)  =  XI  I C  ,  2  )*  ,  2  74  +  X  I  1C.1) 
X|IC.4)=X|1C,4  )- YLDL  *  1  5. 

XI  1C,  5  )=X|  IC,  5  )*  .0  I 


XI  IC,6)=|X|TC,6)*|X{  IC,I)-XMATLI  1  ,1D))  »■  1  XMA  TL 1  I 
1274)*  1SC0.  )/I XI  IC,  1  )♦ XMATLI 2,1 B) *1 .274, 

TSV=ITSI  1C,  IF)  )  4-TSL  )*RNORMI  4) 

RETURN 

2  no  3  1=2,5 

3  XI  1C,  I  )=Q  ,0 


,1  B>  +XMATLI2, IB>* 1 . 


XI  IC,  1  )  =  X|  IC,  1  )♦ XMATLI  1,I0)*RN0RM(|) 

XI  IC,  ft  >=  ;  XI  IC.c.  )*(  XI  IC,  1  )- XMATLI  I  ,  I  B  >  >  *  (  XMA  TL  I  1 
1274)*  150  0. )/(X(  !C,  I)*XMATLI2,IB)*1,274) 

TSV= I TSI  IC,  13  )  4  rSL  )*CNORM|  4) 

RETURN 

ENC 


.1  B)*XMATLI2, IB)*S , 
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Appendix  1  ( Cont 1 d . ) 


Subroutine  (CONVT) 


♦  ONE  MORC  INTEGERS 

SUBROUTINE  C0NVT1TMC1 
DIMENSION  XM( 141 

COMMON  10.  lM.JNIT.JEVNT.JMNIT.MFA.MSTOP.MX.MXC, NCLC  T . NHI ST . 
1N0C.N0RPT.N0T.  NPRM S.NR UN .NR  UNS  « NSTA  T .OUT .SCALE .1  SEED  «T  NOW . 
2TBEG.TFIN.MXX. NPRNT.NCRDR .NEP, VNO<  4 ) 

COMMON  ATR  1  0<  4  1  •  ENO<  4  I  ,  1  NN<  >.  1  .  JCt  LSI  5.22  1  .KRANK14)  .  J  CL  R.  MAX  NO  14)  e 
IMFEC4).MLC141.MLE( 41 .NCELSI 5) .NO  I  41 .PARA Ml  20.4) . QT 1 MEC4 I.SSUMA110, 
251. SUM A( 1 0. S )  .  NAME  I  6  1  .NPRO J.MON  «NDA  V •  N  Y  R 
COMMON  XMATL1J. 61.X1 3.61. XN< 141.XD16I .CRND . I C .IB.NB.TCU131 .TS13.6I 
1. CELT. CU(3 l.OXYl C 1 
MR  ITEfKlPRNT.  101  11  C.  IB 
SCL=  10. 

101  FORMAT (5X.2161 
T=  0  • 

IFC 1H-CI  1.2.2 

1  I F 1 XI  1C. 6 1-2250 . 1 3.  4, 4 

3  IF1X1IC.2114.20.23 

20  1F(X( 1C. 3114.21. 21 

21  JL  =  1 

MR  1TF1NPRNT.  1001T.1X11C.L1.L=1.61 
JN  — 0 

DO  6  1=1.6 

6  XD(  1  1=0. 

DO  15  1=1.14 

15  XM(I !  =  XN<  I  1 
K=2 

DO  A  1=1.2 

XM1JL*  I  1=XM(  JL  +  1  1/1  XI  IC.  I  14X1  IC.214X1  1C  .31  I 
DO  7  J= 1.6 
JL= JL ♦ 1 
JN=JN4 1 

7  XC  1 K  1=  XM 1  JN  !♦  X 1 1C.J14XD1K1 
JL= JL ♦ 1 

JN=JN+ 1 

XO (K 1=XD<K 14XM ( JN I 
6  K= K  4  4 

u=oxvi  ibi 
X  D ( 1 1  =  0. 

X0121=XD(21-.00184U 
XD ( 3 1  =  XDC  2)4 .274 
XD ( 4 1=  XDI 2 14  1- .866 1 
X0(5!=XD( 2141- .191 
XDI61=X0!614.0164U 
DO  9  1=1.6 

9  XI IC.  I1*X1 IC. I  14X01  I 14DELT 
T=T40ELT 
GO  TO  1 

4  TMC=TXSCL 

MR ITE1NPRNT.  100 1 T . 1 XI  I C . L I . L  =  1  .  61 
RETURN 

2  T  =  0  . 

cu< ic»=x« ic, ii4 .a 
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Subroutine  (CONVT)  Cont'd. 


10  IF(»(IC«l)**in2.Uill 

11  u=oxri  IB)  ,  .. 

WR 1T£ ( NPRNT • 100)T. (  X( IC.L)«L-1 «6) 
ARC=.150E-06*U*T**2. 

X(  1C*  l  )  —  X (  IC.  I  >*EXP<-AftG> 

t=t>oelt 
GO  TO  10 

12  TMC*T/SCL 

WR 1TE  < NPRNT.  100)T.(X(  1C.L).L  =  1  .6) 
100  FORMAT < 2X. F8. 3.6» 2X.E » 4. 4> > 

RETURN 

ENC 
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Appendix  2 

Model  Parameters  and  Assumed  Oxygen  Input  Rates 

Coefficients  of  the  Slag  Phase  Model  as  expressed  by  Equation  3.9 


A11 

85 

-.29 

A51 

= 

.085 

A12 

= 

.0005265 

A52 

= 

.000106 

A13 

= 

.008117 

A53 

= 

-.007658 

A14 

= 

.001004 

A54 

= 

-.0004252 

A15 

= 

-.00006333 

A55 

= 

-.00003307 

A16 

= 

.00769 

A56 

= 

-.005419 

C1 

= 

.06282 

C5 

= 

.2167 

k3 

= 

.274 

b5 

= 

.016 

K4 

= 

-.93 

K6 

= 

-.19 

bl 

= 

-.0018 

Assumed  u  =  150  lb/minute 

Coefficients  of  the  Finish  Phase  Model  as  expressed  by  Equation  3.17 

a  =  .15  x  i0"6 
Assumed  u  =  100  lb/minute 


Appendix  3 

Computer  Algorithm  for  the  Minimum  Blow  Duration  Problem 


109 


5TH  CARD  - 


AFTER  CATA  CARD  INSERT  11  =  1  (FOR  360  COMPILER)  OR  1  I  =  2  (FOR  1803  COMPUTER) 

DATA  CARDS  MUST  OE  PUNCHED  AS  FOLLOWS  AND  STACKED  IN  THIS  ORDER:  CO*,,UTER, 

1ST  CARD  -  CU2S  (AMOUNT  ADDCD  FOR  EACH  OF  5  BLOWS!  FORMAT  5F10.0 

2ND  CARD  -  FES  (AMOUNT  ADDED  FOR  EACH  OF  S  BLOWS)  FORMAT  5F19.0 

3RD  CARD  -  FLUX  (AMOUNT  ADOEO  FOR  EACH  OF  5  OLOWS)  FORMAT  5F10.0 

♦th  caro  -  TDi  (initial  value  of  blow  duration),  deltb  (increment  for  blow 

DURATION).  M3  (NUMBER  OF  TIMES  TO  INCREMENT  BLOW  DURATION)  FORMAT 
2F 1 0 . C.  1  1  0. 

5TH  CARD  -  VI  (INITIAL  VALUE  OF  OXYGEN  RATE),  DF.LTV  (INCREMENT  FOR  OXYGEN 
RATEI.M4  (NUMBER  OF  TIMCS  OXYGEN  RATE  IS  INCREMENTED)  FORMAT 
2F  1  0 . 0  •  1  1  0  „ 

D 1  MENS  ION  T(lO.lO.e),  U(10.10.6).X(6),CU2S(5).FES(5>  .FLUX (5) ,XS (6) 

DATA  T.U, X/ 1 206*0 . 0/ 

11=1 

RE A0( I  1. 1 )( CU2S(  1 ) . 1  =  1 ,5) 

READ! 1 1, 1 )(FES( 1 ). 1=1.5) 

READ!  11,  1  1( FLU  X (  1 ). 1  =  1.5) 

1  FOKMAT(5F10.0) 

READ! 1 1.2)TE1, DELTB, MA 
RE  AC (  11,2)  V1.DELTV.M3 

2  FORMAT (2F10.0, 1 10) 

X ( 4  )  =  0 .0 

X(S)  =  1900.00 
X (6 )  =  0 .0 
X(  1)  =  FES(  1  ) 

DO  35  N=  1, 5 
X ( 2 )  =  CU2S(N)  ♦  X ( 2  ) 

X  (  3  )  =  FLUX(N) 

N  1=  1 
N2=  1 

1FIN-1 ) 20 , 20. A 
A  IT  =  T (N1.N2.N ) 

I F ( IT-0131, 31,5 
5  X  (  I  )  =  5  *  N 1 
X (  1  )  =  X (  1  )  -  2.5 
X (  1  )  =  FES ( N  )  ♦  X (  1  ) 

X(S)  =  50 ,+N2 

X(S)  =  X ( 5 )  -  25.0  +  1900.00 
A  =  FES(N)  ♦  CU2SIN)  ♦  FLUX(N) 

B=  X ( 1 )  ♦  X ( 2 )  -FES(N)-CU2S( N) 

C  =  X(  1  )  *  X  (  2  )  ♦  X  (  3  ) 

X (  —  )  =  (X(5)  *  B  ♦  ISO  0.0  *  A)/C 

20  V  =  V  1 

DO  30'M1=1,M3 

TC  =  Tei 

DO  29  M2=  1  ,  M  A 

T0T  =  TB  ♦  TIN1.N2.N1 

CALL  DONOT  (TB.V.X.XS) 

IF (N-5  )  I3CC. 1001,  1C01 
1001  1F(XS(  1)-5.)2|,21.2S 
1000  IF(XS(  1  )-0. 0)29,29,21 

21  1F(XS(5)-22S0.0)22,22,29 

22  00  23  1=1,6 
R  =  5*  1 

IF( XS( 1  )  -R  )  2  A ,  2A,  23 

23  CONTINUE 
2A  S=  1 950  . 

OD  25  J= 1,7 
1FIXSI 5 ) - S ) 26 .  26,25 
25  3  =  S  ♦  50.0 
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Appendi x  3  ( Cont 1 d. ) 


26  NS*N»1 

IT  =  T(  I.J.NSI 
JF  I IT-0 128,  28,  27 

27  IF|T< I, J.NSI-TRTI29,  29.28 

28  TdiJtNSI  =  TBT 
U(  I.J.NSI  =  V 

WR ITE ( 3,  70 ITBT,  V, I  . J.N1.N2 
70  FORMAT (SX.2E12.5.4I6I 

29  TB  =  TE  ♦  OELTB 

30  V  =  V  ♦  OELTV 
IFIN-I |3S. 35. 31 

31  N2  a  N2  ♦  1 

IF (N2-7  14, A, 32 

32  N2  =  1 

N 1  -  N 1  4  | 

I F IN  1- 6  14.4.3S 
35  CONTINUE 
wniTE(3.44l 

441P°'1AT<'  *  BL°'“  T°TAL  BL0“  TIMt  OXYGEN  RATE  ICFES) 

00  50  K=  2, 6 
00  SO  1=1,6 
00  50  J= 1,7 
L  =  K-1 

IT=T(I,J,KI 
IFI  IT-0150. SO. 45 

45  WRITEI3.46IL.TI  l.J.Kl.UH.J.Kl  ,1  ,  J 

46  FORMAT!*  *  .  2X,  I  1 , 5X.F  1  0 . 2, 6X  ,F  1  0.  2.  HX,  12  ,7X,  121 
50  CONTINUE 

CALL  EXIT 
ENC 

SUBROUTINE  OONOTI TEND . U. X, XS I 

DIMENSION  XI6I.XSI6I 

XSI 1 !=X( 1 1- .0O1774L4TEND 

XS  (5I  =  X(  5  M,  015654  L*  TEND 

RETURN 

ENC 


JIT  EM 


Appendix  4 

Computer  Algorithm  for  the  Optimum  Scheduling  Problem 


lDi’JET?iQ«,TTn?;6i:T!l!3,6,,TB<3,6‘*TBM3*6,‘TS0<3*6»*Ta0<3.6I.TD<3 
1  ,  A  1 ,  T  (  3.  6  )  ,  TO  (  .1 , 6  )  ,  A  T(  3,6)  ,LB(  3.6)  .  TO  0(3  .6) 

DATA  TOO/104O.O/ 

"S  *1°  ,NPUT  VALUES  OF  SERVICE  TIME  AND  MINIMUM  BLOW  DURATION 

READ!!,  l  KTSII  I.K).K=1.6> 

RE  AC!  1.  | )  ITQH  |  ,K  I.K=l  .6) 

FORMAT(6F10.0) 

110  CONTINUE 

WR  ITE<  3. 39) 

39  FORMAT  (•  •,  'INPUT:  CONVERTER  riLOW 

IN'  ) 

DO  30  M= 1.3 
00  30  N=  1.6 

MR  ITEI 3. AO )M,N.TSI (M.N I . TR 1 (M.NI 
*0  FORMAT ( •  • , 1 IX , 1 1. 6X, 1 1 , 5X,F 10. 2.5X,FI 0.21 
30  CONTINUE 

ROTATES  INPUT  VALUES  TO  TRY  EVERY  COMBINATION 


SERVICE  TIME  MIN.  BLOW  DURATIO 


OT  T-  0 
00  13 

00 
00 


=  1.6 
13  1  A=  1 , 6 
13  10=1,6 


GEDOQMFL-,|fa  ARRArS  DF  VAL'JES  F0R  EACH  COMOINAT10N  OF  INPUT  DATA 

LC=L+I-1 

GO  TO  < 47, 47, 47, 47,  47.  47 . 4 1 , 42 , 43. 44 ,45 ,46 )  , LC 

41  L  C=  1 

GO  TO  47 

42  LC=? 

GO  TO  47 

43  L  C=  3 

GO  TO  47 

44  LC=4 

GO  TO  47 

45  L C=  5 

GO  TO  47 

46  LC=f 

47  TS<  I.L  )=TSl<  1.  LC  ) 

Ten  ,L  I  =  TB1|  1 ,  LC  ) 

00  62  LA=1,6 
LD=LA  ♦  ia  -  1 

GO  TO  I  54, 54,  54,  54, 54,  54 , 48 , 49 ,50  .51  ,52  ,53  )  ,L0 

48  LD=  1 

GO  TO  54 

49  L  0=  2 

GO  TO  54 

50  L  0=  3 

GO  TO  54 
5  1  L  0=  4 


GO  TO  54 

52  L0=  5 

GO  TO  54 

53  L 0=  6 

54  TS(2,LA)  =  TSH2,L0I 
TB(2,LA)  =  TEIH  2.LD  I 
00  62  LL=  1,6 
LE=LLT IF3-1 

GO  To  16 1 . 6 1 . 6 1 , 6 1  •  6 1 , 61  , 55. 56, 57 ,58. 59 ,60)  ,LE 
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55 

L  £=  1 

GO  TO 

6  I 

56 

LE=2 

GO  TO 

61 

57 

LE=3 

GO  TO 

61 

59 

LE=4 

GO  TO 

61 

59 

LE=5 

GO  TO 

61 

60 

LE=6 

61 

TS  C  3.LL )  =  TSIC  3.LE ) 

62 

T8C3.LLI  =  TBIC3.LF) 

COMPUTES  TOTAL  TIME  FOP  A  COMPLETE  C VC LE  (6  BLOWS  AT  ALL  3  CONVERTERS). 

00  2fl  1M=1,3 
00  26  IN- 1*6 
28  T0( IM,  IN  I  =  0.0 
00  11  K  AN  = 1.3 

00  11  J=1.6 
J 1=J-1 

1  F  ( J  1  -  0  1  3 .  2,  3 
2  J  1  =  6 
1  J2= J ♦ 1 

1FU2-7I5.4.5 
A  J2=I 

5  T( l. J )*TS< 2. J )+TSC 3, Jl+TOC  2, J1  )*TOC  3, J1  I 
TD( 1. J  )  =  T  0  C  l.J  )-TC  I, J) 

IF(TO(  l.J  )  16.6. 7 

6  TOC  1. J  )=0 

7  T(  1. J  )  =  TC  l.J) tTDC  l.J) 

TC2,  J)  =  TM  1.J2IFTSC  3.JM-T0C  1  .  J  ) ♦  TO ( 3  .  Jl  ) 

TCC2.J  )  =  T  BC 2, J )-TC  2.J) 

IFCTOC  2 , J ) )8,8,q 

8  TO  (  2.  J  )=0 

9  T  (2.J  )  =  T(  2.J  I  +Tl)C  2  «  J  ) 

T ( 3 . J I  *T S  C 1.J2I+TSC 2.J2ITTDC  1 , J)  +  TDC  2 , J) 

T0C3.J  )-TBC  3.J  )-  TC  3.J) 

IFCTOC  3, J ) ) 10.  10.  1 1 

10  T0C3. J  1  =  0 

11  TC  3. J l  =  T  C 3. J ) +T0 C  3. J I 
TT  =  0 

00  22  K =1.3 
00  22  L=  1 . 6 

22  TT=TOCK.L)  ♦  TSCK.L)  +  TT 

SFLECTS  MINIMUM  CYCLE  TIME  AND  STORES  ARRAYS  QF  SERVICE  TIME.  MINIMUM  BLOW 
OUR AT  I ON .  ANO  OPTIMUM  BLOW  DURATION 
TT  1= 1  ,/TT 

I  F  C T T I -OTT ) 1 3.  13.  12 

1 2  OTT=TT I 

00  21  M=l,3 
DO  21  N=  1 . 6 
TSOCM,  N )  =  TSCM,N  ) 

TBOCM, N )  =  TBCM,N  I 
TOOCM.N )=TOCM,N  ) 

21  TOCM.N )  =  TCM,N  ) 

10=  I 
I A0=  I  A 

iec=  io 

TT0=1 ./OTT 

13  CONTINUE 
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C  OUTPUTS  RESULTS  OF  OPTIMIZATION.  STATES  OPTIMUM  STARTING  COND:TIONS  FOR 
C  EACH  CONVERTED  AND  fJPTlNUM  TOTAL  CYCLE  TIME. 

W  R  I TE (3*  14)10,  1  AO*  100, T TO 

14  FORMAT ( • 0* • *  A  T  T=0:  CONV.  I  IS  ON  CLOW  •.II.*.  CONV.  2  IS  ON  BLOW 

1  •I'***  C0NV*  3  >S  ON  BLOW  *. !!.•.•/  •  THIS  GIVES  OPTIMUM  TOTAL  TI 
2ME  OF  *,F10.2.*  MINUTES.*) 

C  GENERATES  AN  1N0CX  r0R  EACH  CONVERTER.  To  INDICATE  CORRECT  BLOW  NUMBERS 
00  71  L  =  1 , 6 
L  B (  1  *  L  >  =  L  ♦ 10  -1 
L  B  (  2*  L  )  =  L  ♦  1  AO  -  1 

71  LH I 3.L 1  =  L  ♦  1 00  -  1 
DO  7  7  M=  1  ,  J 

DO  77  N=  1 , 6 
L  =  LOIM.NI 

GO  TO  I  T7,  77,  77,  77,  77,  77,  72.  73 • 74 ,75.76)  .L 

72  L3(M,N 1= 1 
GO  TO  77 

73  LBIM.N  >  =  2 
GO  TO  77 

74  LBIM.N  1  =  3 

GO  TO  77 

75  LBIM.N)  =  4 
GO  TO  77 

76  LBIM.N )  -  5 

77  CONTINUE 

.  OUTPUTS  AN  ARRAY  OF  OPTIMUM  BLOW  DURATIONS 
WR  1TEI 3.  100 ) 

100  FORMAT  1 •  1 *, 1 6X.  'OPTIMUM  OLOW  DURATIONS* //•  *.15X.*BL0W  1  BLOW  2 

1  eLOW  3  BLOW  4  BLOW  5  BLOW  6*) 

DO  103  1=1,3 

WRITE  13.  101  >1.  ITOI  1,K),K~1«6) 

101  FORMAT  I • 0 • , •  CONVERTER  * . I  1 . 61  2 X ,F 7 . 2 ) ) 

103  CONTINUE 

GENERATES  A  TIME  SCHEDULE  FOR  OPTIMAL  CRANE  SERVICE  FOR  CONVERTERS. 

DO  B5  1=1.3 
DO  BE  J  =  1 , 6 

85  AT  I  I , J  )  =  0.0 
DO  82  KEN=  1 , 3 
1 FIKEN- 1  196. 96,98 

96  AT  I  1.  1  >=0.0 

AT  12,  /  |  =  T SO  I  1,1) 

AT  I  3,  1  )  =  TSOI  1 .  1  )  ♦  TS0I2.1) 

DO  97  1=  1,  3 
DO  9  7  J  =  2 , 6 

97  ATII.J)  =  ATIl.J-l)  ♦  TOIl.J-i)  ♦  TSOII.J-l) 

GO  TO  99 

98  CONTINUE 

DO  70  K—  1,3 
DO  70  L  =  1  •  6 
L3=L-1 

I F  IL  3-0 170.  95.  70 
95  L  3  =  6 

70  ATIK.L)  =  ATIK.L3)  ♦  TOI<C,L3)  ♦  TS0IK.L3) 

99  CONTINUE 

OUTPUTS  CRANE  TIME  SCHEDULE. 

WR  1TEI 3. 80 ) 

"^O^ER^cTr,^  OLOW  DUR 

DO  82  J=  1 , 6 
DO  82  1=1,3 

WRITE  I  3* 81  1ATI  I , J ),  1 ,LD<  I . J) ,TO<  1  ,J) , TSOI  1  II 

. 

CALL  FX1T 
ENC 
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Appendix  5 

Computer  Algorithm  for  the 
Converter  Furnace  Control  Problem 


•CNu  WOHD  INTEGERS 
♦LIST  SOURCE  PROGRAM 

SUEFQiJT  INC  TONUO(  TCND.U.  XX,  XTfc  ST.  T) 

DIMENS  ION  XMI|4>,A(H»,x*(6).XS(6|,XD<6»,0UMYXt6J 

1  -0.6333E-04,  O.7t'ICE-O2°02A?3Pe  ’r  ;5n65E-33 ,0*BI  1  7F‘02,°  •,004E-02  , 

**  blo“  r°°  ^ 

•  XM(|>  is  Thfc  A  I  I  , J )  COEFFICIENT' 

OELT=5.0 


c., 
c .  , 


U  =  25  . 

•  •..DUMMY  ARRAY  TO  PRESERV  All 
DO  125  NO=  |  , A 

a<nci=xm<nci  i 

125  CONTINUE 

1 1  T=a. 


FOR  CALCULATING  All 


PRIME 


DO  124  N’=|,6 
XS ( NP )=  X x ( NP  I 
124  CONTINUE 
13  JL  =  0 


CO  19  1L3 1 • 6 
XC(ILI=C. 

19  CONTINUE 
T=  T  +  CfcL  T 
CO  Eg  IY=1,6 
dumyxi  iyi=xs( iy> 

69  CONTINUE 
K=  I 

DO  16  1-1,2 

XM  <  JL  «■  I)  =  A  <  JL  +  I  >/(  X  S (  1)  +  XS(2)+X'(3)I 
C3  I  5  J= | , 6 
JL=JL+ I 

XD(X  )=XM(JL  )*XS( JI  +  XDIK I 

15  CONTINUE 
!L=JL* I 

XD(K  )=XD(X  )«-XM(  jl  ) 

K  =  K44 

16  CONTINUE 

XD(  1  »=XC(  1 »-  .GC  1F*L 
XD(3)=XC( I }* .274 
xo(4)=xci i i»(- .gtA) 

XO  (  5  1=  XD  (  5  >♦  .0,16*0 
XD(S)=X0(1)*(-,19) 

DO  25  1=1,6 

*S<  I  >=  XS<  1 l+XDI  1  >*OLLT 

25  CONTINUE- 

I F ( X  S (  l»-XTEST)77,  77,|CC 
133  IF(T-TENCI13,7g,7g 
T9  U=  U  ♦  5  . 

GO  TO  II 

77  1F(T-(TcND-5.)»44.  ICl,  1C1 

10  1  I F ( T - ( TEND *5 . I  I  1 C2 ,  1 02 , 76 
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. LINEAR  INTERPOLATION 

M  CO  61  I  Z“  I i 6 

XS(  I  7  >  =  <  <S(  I  7 ) 40  00  VX<  17)  1/7. 

61  C  JNT INUF 

r-T-( OtLT/? . ) 

102  *h|TF(’,2|T,u,(XS(  l),r  i,6) 

2  FORM  AT (HE1  1  .4 > 

he  ruRM 

44  W  f.  I  T  E  (  3  .  4  6  ) 

46  FORMAT  (•  WHAT  H.\nP|.Nfr?  M|SStl>  P)  CULT.*) 
(■iTIIRM 
•NC 
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